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ABSTRACT 

The first three months of sky-survey operation with the Fermi Gamma Ray Space 
Telescope (Fermi) Large Area Telescope (LAT) reveals 132 bright sources at jb| >10° 
with test statistic greater than 100 (corresponding to about lOcr). Two methods, based 
on the CGRaBS, CRATES and BZCat catalogs, indicate high-confidence associations 
of 106 of these sources with known AGNs. This sample is referred to as the LAT Bright 
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AGN Sample (LB AS). It contains two radio galaxies, namely Centaurus A and NGC 
1275, and 104 blazars consisting of 57 flat spectrum radio quasars (FSRQs), 42 BL Lac 
objects, and 5 blazars with uncertain classification. Four new blazars were discovered on 
the basis of the LAT detections. Remarkably, the LBAS includes 10 high-energy peaked 
BL Lacs (HBLs), sources which were so far hard to detect in the GeV range. Another 
10 lower-confidence associations are found. Only thirty three of the sources, plus two at 
|b[ <10°, were previously detected with EGRET, probably due to the variable nature 
of these sources. The analysis of the gamma-ray properties of the LBAS sources reveals 
that the average GeV spectra of BL Lac objects are significantly harder than the spectra 
of FSRQs. No significant correlation between radio and peak gamma-ray fluxes is 
observed. Blazar log N - log S and luminosity functions are constructed to investigate 
the evolution of the different blazar classes, with positive evolution indicated for FSRQs 
but none for BLLacs. The contribution of LAT-blazars to the total extragalactic 7-ray 
intensity is estimated. 

37 Subject headings: gamma rays: observations — galaxies: active — galaxies: jets — BL 

38 Lacertae objects: general 



Introduction 



40 The Gamma ray Large Area Space Telescope {GLAST ) was launched on 11 June 2008, and re- 

41 named the Fermi Gamma Ray Space Telescope shortly after entering its scientific operating mission, 

42 which began on 11 August, 2008. The Large Area Telescope (LAT) on Fermi provides an increase 

43 in sensitivity by more than an order-of-magnitude over its predecessor EGRE T, the Energetic 



44 Gamma Ray Experiment Telescope on the Gompton Gamma Ray Observatory (iThompson et al 



I993I ) , and the Italian S pace Agency Satellite AGILE (Astro-rivelatore Gamma a Immagini Leg- 



46 gero; iTavani et al.ll2008l ). In sky survey mode, the LAT observes all parts of the sky every 3 hours, 

47 providing effectively uniform exposure on longer timescales. 

48 One of the major scientific goals of the Fermi Gamma Ray Space Telescope is to provide new 

49 data about 7-ray activity of AGNs. Rapidly varying fluxes and large luminosities of extragalactic 

50 7-ray sources are best explain ed if the 7 rays are emitted from collimated je ts of charged particles 

51 moving at relativistic speeds (jBlandford &: Reeall978l : iMaraschi et al.lll992l ). Fermi-LAT observa- 

52 tions will help determine how these particles are accelerated, where the gamma rays are emitted, 

53 what the energy and power budgets of the supermassive black-hole engines are, what this says for 

54 the fueling and growth of black holes, and the reasons for the differences between radio-loud and 

55 radio-quiet AGNs, and FSRQs and BL Lac objects. These are jus t a few of the ques tions that 

56 7-ray AGN studies with the Fermi-LAT are helping to answer (see lAtwood et al.l |2009| . for more 

57 discussion of these goals). 



In a companion publication (jAbdo et al.l l2009d ). 132 bright sources at |6| > 10° with test 



- 5 - 



59 statistic (TS) > 100 are found in the preliminary three month Fermi all-sky survey. As expected 

60 from the EGRET legacy, a large fraction of these sources are AGNs. Detailed results of the subset 

61 of the Fermi bright source list that are associated with AGNs are presented here. 

62 Sixty-six high-confiden ce blazars are listed in the Third EGRET catalog of high-energy gamma- 

63 ray sources (3EG catalog; iHartman et al.Ml999l ). with the majority of them, ~ 77%, identified 

64 as flat spectrum radio quasars (FSRQs), and the remaining ~ 23% identified as belonging to 

65 the BL Lac c l assic The recently released catalog of high-confidence AGILE gamma-ray source^ 

66 ( Pittori et al. 20081 ) shows a somewhat higher percentage of BL Lacs. Unlike AGN surveys at 



67 optical or X-ray energ ies, in which the majority of AGNs are radio quiet (e.g., Idella Geca et al 



19941 : llvezic et al.ll2002l ). all AGNs detected at > 100 MeV energies are also significant radio sources. 



69 This includes the 3EG and AGILE blazars, which are so far identified with flat spectrum (radio 

70 spectral index Or > —0.5 at GHz frequencies) ra dio-loud AGNs, and most show superluminal 



71 motion Oorstad et al. 



2001 



Kellermann et al.ll2004l ). Moreover, the redshift distribution is broad, 

72 with the largest redshift AGN known in the 3EG catalog at z = 2.286. 

73 Here we present a source list of bright AGNs found in the set of the 132 bright LAT sources at 

74 |b| >10°. Identiflcation of variable 7-ray sources with blazars depends on the statistical likelihood 

75 of positional a ssociation and correlated var iability of the 7-ray emissions with lower-frequency 

76 radiations (e.g. ISowards-Emmerd et al.ll2003l ). The 106 sources having high-confldence associations 

77 with known blazars and radio-galaxies constitute the LAT Bright AGN Sample (LBAS). Included in 

78 this list are mean fluxes, weekly peak fluxes, spectral indices, locations, and variability information. 

79 Only sources with confldence levels greater than lOcr are retained in the LBAS. This list is not, 
BO however, complete, as we already know of many more sources at lower significance. The limiting 
Bi flux depends on both the source sky location and the spectral hardness. 

B2 In Section 2, observations with the LAT, analysis methods, and the source detection procedure 
B3 are presented. Section 3 describes the association method and gives the list of bright Fermi-hAT 
B4 detected blazars. Key properties of the LBAS, including flux and spectral index, are presented 

85 in Section 4. The LBAS is compared with EGRET blazars in Section 5. Section 6 considers 

86 the radio/gamma-ray connection. Population studies, including source types and redshifts, are 

87 presented in Section 7, where the log N - log S flux distributions and luminosity functions of the 

88 LBAS are constructed. The results are discussed in Section 8, including implications of the results 

89 for blazar evolution. We summarize in Section 9. 

90 In the following w e use a ACDM cosmology with values given within la of the WMAP results 

91 (jKomatsu et al.l l2008l ) . namely h = 0.71, = 0.27 and = 0.73. Here the Hubble constant 



^In contrast to the prominent optical emission lines found in FSRQs, BL Lac objects are radio-loud, rapidly 
variable sources displa ying nearly featureless continua with emission-line equivalent widths < 5 A (for review, see 
Urrv fc PadovanilllQQSh . 



^http://www.asdc. asi.it/agilebrightcat/ 
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Hq = lOO/i km s ^ Mpc ^ is used. 



2. Observations with the Large Area Telescope 



94 The Fermi-LAT is a pair-conversion gamma-ray telescope sensitive to photon energies greater 

95 than 20 MeV. It is made of a tracker (composed of two sections, front and back, with different 

96 capabihties), a calorimeter and an anticoincidence system to reject the charged-particle background. 

97 The LAT has a large peak effective area (~ 8000 cm^ for 1 GeV photons in the event class considered 

98 here), viewing ~ 2.4 sr of the full sky with excellent angular resolution (68% containment radius 

99 ~ 1° at = 1 GeV for the front section of the tracker and about a factor of 2 larger for the back 

100 section). A fu l l desc ription of the LAT instrument and its predicted performance are reported in 



Atwood et al.l (j2009l ). During the first year, the telescope operates in sky-survey mode observing 



102 the whole sky every 3 hours. The overall coverage of the sky is fairly uniform, with variations of 

103 around ~15% around the mean value. 

104 The LAT data used here were collected during the first 3-month all-sk y survey, from August 4 

105 to October 30 2008. We refer to the companion paper (lAbdo et al.ll2009a ) for a full description of 

106 the data selection and analysis. In order to avoid background contamination from the bright Earth 

107 limb, time intervals where the Earth entered the LAT Field-of-View (FoV) were excluded from this 

108 study (corresponding to a rocking angle < 47 deg). In addition, events that were reconstructed 

109 within 8° of the Earth limb were excluded from the analysis (corresponding to a zenith angle cut 
no of 105°). Due to uncertainties in the current calibration, only photons belonging to the "Diffuse" 

111 class with energies above 100 MeV were retained. These photons provide the purest gamma-ray 

112 dataset. The energy range was even more restricted in the source detection and spectral fitting 

113 analyses described below, where only photons with E >200 MeV were selected. The list of sources 

114 reported in Tables 1 and 2 was obtained as the r esult of the source detection, localization and 

115 significance estimate analyses described in detail in lAbdo et al.l (j2009d ). 



116^ ^The source detecti on step made use o f two wavelet algorithms, {mr_filter) (jStarck &: Pierre 



19981 ) and {PGWAVE) (jCiprini et al.l 120071 ). The algorithms were run independently for different 



118 energy bands associated with different localization power and the results were cross-checked. The 

119 positions of the sources for which the detection significanc e was above thresho ld (4(t) were then 

120 refined using {pointfit), a simplified likelihood method (see lAbdo et al.l (|2009d )). This algorithm 

121 uses photons with E>500 MeV and returns the optimize d sky position as we ll as an estimate of the 

122 error radius for most detected sources. As discussed in lAbdo et al.l (|2009d ). the final error in the 

123 source position was estimated by multiplying the error radius returned by the algorithm by a factor 

124 close to 1.4 and adding 0.04° in quadrature (estimated from the residuals between the estimated 

125 and expected position of Vela). The 95% confidence error radius was then evaluated assuming a 

126 2-D normal distribution. 



127 



To better estimate the source significance, we used the maximum likelihood algorithm imple- 
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128 mented in {gtlike) a tool that is part of the standard Fermi-h AT ScienceTools software packag^. 

129 The flux, photon index and test statistic (TS) of each source in the energy range 0.2-100 GeV were 

130 determined by analyzing regions of interest (ROI) typically 15° in radius. The model of the ROI 

131 used to fit the data was built taking into account of all the sources detected within a given ROI. 

132 The isotrop i c back ground and Galactic Diffuse background models used in the fit are discussed in 



i33 lAbdo et al.l ( 2009d ). Each source was modeled with a simple power law (kE""^) for photons E> 200 

134 MeV. The flux [E>100 MeV] (Fioo), which is conventionally reported, was then calculated with the 

135 fitted parameters. This flux will be used throughout this paper. The spectral energy distributions 

136 of some bright sources show clear evidence for a break or curvature. A fit with a single power law 

137 function is certainly not the most appropriate choice for these sources but the resulting photon 

138 index does reflect the spectral hardness. A more detailed spectral analysis of the LBAS sources 

139 is beyond the scope of this paper. The source fluxes were also estimated by fitting independent 

140 power law functions in two energy bands (0.1-1 GeV) and (1-100 GeV) and summing up the two 

141 obtained fluxes. These fluxes (Eg 5 in Table 3) are the same as those reported in the Fermi bright 



142 source list paper (lAbdo et al.ll2009d ). Eor most sources, the fluxes obtained by the two methods 



143 are consistent within 30%. 



144 The same procedure was applied to generate weekly light curves (spanning a 12-week period). 

145 Erom those, the weekly peak flux as well as a variability index (corresponding to a simple 

146 criterion) were derived. The variability tag reported in this paper is set for sources associated with 

147 a probability of being constant lower than 1%. A few representative light curves are displayed in 

148 Eig. [H 

149 This analysis was performed with the preflight instrument response functions (P6_V1). In 

150 flight, the presence of pile-up signals in the LAT tracker and calorimeter left by earlier particles 

151 was revealed in periodic-trigger events. This feature leads to a reduction of the real acceptance as 

152 compared to the predicted one as fewer events pass the rejection cuts, most notably for low-energy 

153 photons. The magnitude of this reduction is still under investigation, but the fluxes reported here 

154 may be lower than the true ones by as much as 30% and the photon indices greater than the true 

155 ones by as much as 0.1 (true spectra could be softer by 0.1 unit in the photon index). Because 

156 of the current uncertainty, no correction has been applied to the results. This uncert ainty applies 



157 unifor mly to all sources. Our relative errors are much smaller (about 3% on the flux. lAbdo et al 



2009d ) . With the accepta nce used in this a nalysis, the measured fluxes of the 3 bright pulsars. 



159 Vela, Geminga and Crab (jAbdo et al.l l2009d ) are found to be compatible within 11% with those 

160 reported in the 3EG catalog. 

161 Eig. [2] shows the 3-month flux sensitivity for TS=100 and a photon index=2.2 as a function of 

162 the sky position, calculated by a semi-analytical, maximum likelihood estimate of the significance. 

163 This estimate takes the actual exposure, the PSE and the different backgrounds (galactic diffuse, 

164 extragalactic diffuse and instrumental) into account. The limiting flux is higher at low galactic 



http://fermi.gsfc. nasa.gov/ssc/data/analysis/documentation/Cicerone/ 
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165 latitude due to a higher galactic diffuse background and close to the celestial south pole (1 

166 b ~ —27°) where the exposure is lower. 



302°, 



167 The final result of the detection analysis is a list of 205 sources with a (T S > 100, ^ lOcr), 

168 composing the LAT Bright Source (OFGL) list (see Table 6. in lAbdo et alJ[2009cl ) . For comparison 

169 31 so urces detec ted by EGRET have a sign ificance greater than lOo" in the 3EG (jHartman et al 



19991 ) and EGR (jCasandiian Greniei]|2008l ) catalogs. Of these, only 13 were detected at |6| >10°. 



171 In the OFGL, a total of 132 sources, including 7 pulsars, are present at \b\ >10°. We have explored 

172 the possibility of associating AGNs with the 125 remaining sources. 



Source association 



174 Any source association procedure primarily relies on spatial coincidence. Fig. [3] shows the 

175 95% error radius vs (TS) for the sources considered here. This radius depends on both the flux 

176 and the photon index, with a mean of 0.14°. For comparison, the average corresponding radius 

177 for the blazars in the 18 month EGRET sky survey is 0.62°. Of the 186 |6| > 10° 3EG sources, 

178 66 (35%) had "high" (but unspecified) confidence positional associations with blazars in the 3EG 

179 catalog. Another 27 positional coincidences were noted at l ower significance. Although subsequent 



work (e.g. iMattox et al. 



2001 



Sowards-Emmerd et al.ll2003l ) did find additional associations, ~40% 



131 of the high-latitude 3EG sources remained unidentified. 

132 Although the LAT localization accuracy is much better than those of previous gamma-ray tele- 

183 scopes, it is not good enough to enable a firm identification of a LAT source based solely on spatial 

184 coincidence. For the LAT, a firm identification is assumed only if correlated variability is observed 

185 at different wavelengths. In order to find associations between LAT sources and AGNs, two differ- 

186 ent approaches were pursu e d. Th e first method is based on a procedure similar to that developed 

187 bv lSowards-Emmerd et al.l (|2003l ) for associating EGRET blazars with radio counterparts using an 

188 observational figure of merit (FoM). The seco nd one is based on the c alculation of source associatio n 
139 probabilities following a Bayesian approach (jde Ruiter et al.l 119771 : [Sutherland Sz SaundersI Il992l ) , 
190 similar to that used by Mattox et al.l (l200lh to associate EGRET sources with radio sources. This 



191 method is described in Abdoetal 



(2009d) 



192 Several catalogs were used by the two association methods, the most important on es being 

193 the Combined Radio All- Sky Targeted Eight GHz Survey f CRATES: bealev et al1l2007l ^ catalog 

194 and the Roma-BZCA'^ ( Massaro et al. 2007). The CRATES catalog contains precise positions, 

195 8.4 GHz flux densities, and radio spectral indices for more than 11,000 flat-spectrum sources over 

196 the entire |5| > 10° sky. The Roma-BZGAT is a master list of blazars based on an accurate 

197 examination of literature data and presently includes about 2700 sources, all observed at radio 
193 and optical frequencies and showing proper characteristics of blazars. Sources are classified as BL 



*http: / /www. asdc.asi.it /bzcat 
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199 Lacertae objects (BZB), flat spectrum radio quasars (BZQ) or as blazars of uncertain type (BZU). 



3.1. The Figure-of-Merit Method 



201 The figure of merit (FoM) approach requires a large, uniform all-sky sample of radio sources 

202 from which to draw; for this purpo se, we use the Combined Radio All-Sky Targeted Eight GHz Sur- 

203 vey (CRATES: iHealey et alJl2007l ) catalog. In order to quantify the correlation between CRATES 

204 sources and LAT detections, we compare the average number of positional coincidences between 

205 LAT sources and CRATES sources to the number of positional coincidences between LAT sources 

206 and sources drawn from 1,000 randomized simulations of the radio sky. We count as a positional 

207 coincidence any occurrence of a radio source (real or simulated) within twice the 95% error radius 

208 of a LAT source, and we generate the simulated radio skies by scrambling the Galactic coordinates 

209 of the CRATES sources while keeping their radio flux densities, spectral indices, and counterpart 

210 RASS fluxes intact. 

211 We define the excess fractional source density of radio/7-ray matches as n = l—(iVi.and /-^crates) 

212 and we compute this quantity in bins of radio fiux density 5* 8.4 at 8.4 GHz, rad io spectral index a, 

213 and X-ray flux Fx from the ROSAT All-Sky Survey (RASS: l\^es et al.lE999l ^. These functions— 

214 n{S), n{a), and n{Fx) — constitute the counterpart spectral energy distribution (SED) components 

215 of the FoM. The flnal component is the dependence on the offset between the radio position and 

216 the LAT position, which we model simply as n 



pos 



1 — CL, where CL is the confldence limit 
217 of the LAT localization contour passing through the radio position. The FoM is then given by 
213 100 X n{S) X n{a) x n{Fx) x npos. To evaluate the signiflcance of the FoM, we again generated, in 

219 the manner described above, 1,000 random simulations of the radio sky and computed the average 

220 distribution of FoM. We compared this to the distribution of FoM for the real CRATES sky by 

221 again computing the excess fractional source density as a function of FoM. This fractional excess 

222 can be directly interpreted as a probability Pi of radio/7-ray association for source i, giving an 

223 immediate mapping from FoM to association probability for each individual source (i.e., 1 — Pj is 

224 the probability of a false positive association). We flnd that 1,000 simulated skies result in sufficient 

225 statistics in each FoM bin to ensure that the mapping is robust. Very similar results are obtained 

226 with 10,000 simulations. 

227 The results of this association procedure are shown in Table 1 and Table 2. Most of the 



228 assoc iated radio sources are in the Candidate Gamma-Ray Blazar Survey (CGRaBS; iHealey et al 



20081 1. an optical survey of the 1,625 CRATES sources that were most similar in their radio and 



230 X-ray properties to the 3EG blazars. Optical spectroscopy of the sources with unknown redshifts is 

231 ongoing. We also considered the possibility of an association with a non-CRATES radio source when 

232 no CRATES association was found. Indeed, a FoM can be computed for any object for which the 

233 necessary radio data are available. Thus, for those LAT sources without CRA TES associations, we 

234 drew candidate counterparts from the 1.4 GHz NRAO VLA Sky Surve y fNySSMCondon et al.lll998l ) 

235 or the 843 MHz Sydney University Molonglo Sky Survey (SUMSS: Mauch et aLlbood ). searched 
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236 NED for archival 8.4 GHz data, and calculated the FoM for each candidate. These procedures 

237 find high-confidence {P > 0.90) associations for 101 of the 125 non-pulsar sources in the OFGL 

238 list with \b\ > 10° for an association rate of 81%. We also find low-confidence FoM associations 

239 (0.40 < P < 0.90) for 14 more sources, bringing the total association rate to 92%. Thus, the 

240 radio-bright blazar population continues to dominate the extragalactic sky. 

241 The individual association probabilities can be used to estimate the number of false positives 

242 in a given sample: if the probabilities Pi are sorted from highest to lowest, then the number of false 

243 positives in a sample of k sources is A^'faise ~ J2i=i ~ ^i)- Among the high-confidence associations, 

244 there are ~3 false positives, and less than one of the 74 most probable associations should be false. 

245 We also studied the power of the FoM analysis to reject a blazar association for a LAT source. 

246 We considered NVSS/SUMSS sources in the direction of the unassociated LAT sources and com- 

247 puted the FoM that each source would have if (A) it were as bright as the 4.85 GHz flux density 



243 upper limit fro m the Green Bank 6 cm s urvey (GB6: lGregorv et al.lll996l ) or the Parkes-MIT-NRAO 



249 survey (PMN; [Griffith &: Wrightlll993l ) (unless the source had an actual GB6/PMN detection, in 

250 which case we used the measured flux density) or (B) its radio spectrum were as severely inverted 

251 as a = -1-0.75 between 1.4 GHz and 4.85 GHz, whichever constraint was tighter. From the low- 

252 frequency radio spectrum (or upper limits), we extrapolated the implied 8.4 GHz flux density. If 

253 the resulting FoM indicated that the source could conceivably be a flat-spectrum blazar, then we 

254 drew no conclusion, but if we found that the "best-case" association probability were 0%, then 

255 we concluded that the LAT source was not associated with any typical member of the population 

256 of flat-spectrum blazars, and we refer to such cases as "anti-associations." Note that the spectral 

257 index a = -1-0.75 is an extremely conservative cutoff. The most inverted radio spectrum for any 

258 actual association has a < 0.65. We are able to secure anti-associations for 10 sources. In fact, five 

259 of these turn out to be high-latitude LAT pulsars and pulsar candidates. This shows that, given a 

260 reliable LAT error circle, the FoM analysis is capable of indicating definitively that a source is not 

261 a blazar 



262 3.2. Summary of association results 

263 The combination of the FoM (described above) and positional association methods yields 

264 a number of 106 high-confidence (P > 0.90) associations (constituting the LBAS) and 11 low- 

265 confidence (0.40< P <0.90) associations listed in Table 1 and 2 respectively. Simple extrapolation 

266 of these numbers implies that the LAT should be detecting some 20-25 blazars through the Galac- 

267 tic plane at \b\ < 10°. Indeed, several have already been located, e.g. OFGL J0036. 7-1-5951 (lES 

268 0033-^595), OFGL J0730.4-1142 (PKS 0727-11), OFGL J0826.0-2228 (PKS 0823-223), OFGL 

269 J1802. 6-3939 (PMN J1802-3940), OFGL J1833.4-2106 (PKS 1830-211). A more complete search 

270 for Galactic background blazars, incorporating spectrum, variability and multiwavelength proper- 

271 ties is in progress. 
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272 Tables 1 and 2 report, for each source, the LAT name, the name of the associated source 

273 based on the FoM method, the value of the FoM parameter and its probability, the name of the 

274 positionally associated source and its probability, the redshift and the AGN class. Fig. U] shows 

275 the sky location of the LBAS AGNs. 

276 One source, OFGL J10340+6051 reported in Table 1, merits special comment. Two radio 

277 associations were found by the FoM method for this 7-ray source, one with very high probability and 

278 one with lower, but still significant, probability reported in Table 2. Although the high-probability 

279 source likely dominates the 7-ray emission, it is entirely plausible that the low-probability source 
2B0 contributes non-negligibly to the total 7-ray flux. We believe that as the LAT detects more sources 
2B1 and confusion of the 7-ray sky increases, the power of the FoM formalism will become increasingly 
2B2 important to the identification of multiple lower-energy counterparts of complex 7-ray sources. 

2B3 Fig. [5] shows the overall, normalized angular separation distributions for both sets of sources 
2B4 (i.e. high- and low-confidence associations). The solid curve corresponds to the expected distribu- 
2B5 tion (x^ distribution with 2 d.o.f. ) for real associations, the dashed one for accidental associations. 

286 This figure provides confidence that most associations are real. From this figure, it appears that the 

287 1.4 correction factor applied to the error radius is somewhat overestimated. This overly conservative 

288 factor will be significantly reduced with additional analysis updates. 

289 Four new blazars were discovered. Two of these, CRATES J1012+2439, and CRATES J1032+6051 

290 were classified as FSRQ blazars while CRATES J0144+2705 is a BL Lac. The classification of these 

291 three sources was made on the basis of the broad lines observed in their optical spectrum obtained 

292 after the LAT detection (Shaw et al., in preparation). The forth new LAT detected blazar is CLASS 

293 J1054+2210. Its classification as a BL Lac object was made possible by the analysis of its optical 

294 spectrum available at the SDSS on-line archive. As discussed above, CRATES J1032+6051 is the 

295 source which has a low probability to be associated with OFGL J10340+6051. 

296 The other sources listed in Table 1 and 2 were classified as FSRQ or BLLac following the 

297 i?0TOa-BZCat and CRATES/CGRaBS catalogs. Some sources, which cannot be properly classified 

298 because of the scarcity of available data or which show optical spectra intermediate between those 

299 of BL Lacs, FSRQs or radio galaxies, were assigned to the "uncertain class" ( "Unc" label in the 

300 tables) . 

301 Based on this classification, the LBAS comprises 57 FSRQs, 42 BL Lac objects, 5 blazars of 

302 uncertain type, and 2 radio-galaxies (RGs). The relevant EGRET sample of referenc e corresponds 



303 to that of the 18 mo nth EGRET all-sky survey during Phase 1 of the CGRO mission (IFichtel et al 



304 I1994I : iDermeri I2OO7I ) . This survey had relatively uniform exposure, and contained 60 sources, 46 

305 FSRQs, 14 BL Lacs. BL Lacs make up 40% of the LBAS blazars, a fraction significantly higher than 

306 found with EGRET (23%). The detection of hard sources (BL Lac objects, see below) by the LAT 

307 is intrinsically favored over soft ones (FSRQs). This is partly due to the strongly energy-dependent 

308 PSF. The larger bandpass and higher energy for the peak sensitivity (in the ~ 1-5 GeV range) of 

309 the LAT as compared to EGRET adds to this effect. 
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310 Eleven LBAS sources are associated with blazars already detected in the TeV energy range by 

311 the ground based imaging air Cherenkov telescopes. Among these, 7 are classified as high-frequency 

312 peaked BL Lacs (HBLs): lES 1011+496, Mrk 421, PG 1553+11, Mrk 501, lES 1959+650, PKS 

313 2005-489 and PKS 2155-304; 3 are low-frequency peaked BL Lacs (LBLs): 3C 66A, W Com 

314 and BL Lac and one is a FSRQ: 3C 279. These 11 sources represent more than 50% of the 

315 TeV blazars detected so far (21). The results of simultaneous observations that cover the optical, 

316 X-ray, and high energy gamma-ray bands (LAT and H.E.S.S.) of PKS 2155-304 are reported in 



Aharonian et al.l (j2009l ). Another three HBLs in the LBAS are not yet detected in the TeV range: 



318 KUV00311-1938, 1ES0502+675, B3 0133+388. A total of 10 HBLs are thus present in the LBAS, 

319 a remarkable feature given that sources in this class were difficult to detect in the GeV range. Many 

320 of these sources were not particularly flaring at other wavelengths during the period of observation. 

321 We compared the broad-band (radio, optical. X-ray) properties of our sample of Fermi-LAT 

322 detected blazars with those of the known blazars listed in the i?oma-BZCat catalog and found that 

323 the broadband properties of the Fermi-LAT detected BL Lacs and FSRQs are consistent with the 

324 parent population of FSRQs and BL Lacs. This is illustrated in Fig. [6] displaying the soft X-ray flux 

325 vs radio flux density (at 1.4 GHz) diagram for the Fermi-LAT blazars and the full blazar catalog. 

326 The LBAS includes 13 sources (10 FSRQs and 3 BL Lacs) that were detected in a flaring state 

327 promptly announced to the community through Astronomical Telegrams. Among these, OFGL 

328 J2254. 0+1609, associated with 3C 454.3, is the brightest gam ma-ray extra-galact ic source observed 

329 in the 3-month Fermi-LAT survey and is studied in detail in lAbdo et al.l (j2009al ) . 



330 The Fermi-LAT has discovered gamma-ray emission from a source having an high-confidence 

331 association with NGC 1275, the supergiant elliptical galaxy at the center of the Perseus galaxy 

332 cluster. EGRET observations yielded only an upper limit to the NGC 1275 gamma-ray emission . 

333 All the details about the gamma-ray properties of this source will be reported in lAbdo et al.l (j2009bl ). 



334 Cen A is the nearest rad io galaxy to us and it w as one of the few radio galaxies associated with 

335 a 3EG source (J1324— 4314; ISreekumar et al.lll999l ). It is included in the LBAS and the position 

336 of its nucleus is well inside the 95% confidence error radius of the source OFGL J1310.6— 4301. 



337 The measured Fermi flux is Finn ^ 2.3 x 10" '^ ph cm s , about a factor of 2 greater than that 



338 measured by EGRET (jSreekumar et al.lll999l ) 



339 Recently, two more sources reported in the 3EG catalog were tent atively associated with 

340 radio galaxies , 3C 111 (Hartman et al. 2008), and possibly NGC 6251 (jlVlukherjee et al.l bood : 
34i lFoschini et al.ll2005l ). These objects are not LBAS sources but the number of radio-galaxies de- 



342 tected at high-energy is expected to increase in the near future as more data accumulate. 

343 Table 4 lists the 33 sources associated with 3EG sources (two more located at |6| <10° were 

344 also incorporated). Three bright EGRET blazars associated with 0827+243, PKS 1622—297 and 

345 1730-130 (NRAO 530), whose average EGRET fluxes are in the range of (25 - 47) xlO"^ ph(S > 

346 100 MeV) cm^^ s^^ do not appear in the LBAS. Presumably, these blazars are simply in a lower 

347 flux state than when EGRET was in operation. These 3 sources are also among the 22 sources in 
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343 the pre-launch LAT monitored lislo. Of these 22 sources, 17 have high-significance LAT detections 

349 in the first 3-months of data. The remaining two monitored sources (H 1426+428, lES 2344+514) 

350 did not have previous 3EG detections and thus were not expected to be very bright GeV sources. 

351 We note that the LBAS object B2 0218+35 is a well-known gravitational lens. The source 

352 PMN J0948+0022, associated with OFGL J0948.3+0019, has a flat radio spectrum, but shows an 

353 optical spectrum with only narrow emission lines, making it an "uncertain" -type object in the 

354 Roma-BZCat. 

355 4. Gamma-ray properties of the LBAS 

356 4.1. Introduction 

357 Table 3 lists the key properties of the 116 sources associated with AGNs (sources with low- 

358 confidence associations are in italics): the name, equatorial and galactic coordinates, the TS pa- 

359 rameter measuring the significance of the detection, the photon index (F), the photon flux Fioo, 

360 the weekly peak flux, the photon flux F25 and the variability flag. The uncertainties are statistical 

361 only. From Table 3 (last column), 40 FSRQs (70%), 12 BL Lacs (29%) and 1 Uncertain blazar 

362 (OFGL J0714. 2+1934) present in the LBAS show evidence for variability. The observed variability 

363 for FSRQs is thus higher than for BL Lacs. One must be careful in interpreting this result as the 

364 flux distributions are different for the two classes (see Fig. [7]), making the detection of variability 

365 easier for FSRQs. An in-depth variability analysis of the LBAS is beyond the scope of this paper. 

366 Table 4 gives similar parameters for the subset of 35 sources (including both high-confidence 

367 and low-confidence associations, plus two at |b| < 10°) corresponding to 3EG sources. This subset 

368 will be discussed in more detail in section 5. 

369 The source photon index is plotted as a function of the flux in Fig. [71 It is already visible 

370 in this figure that the photon indices of BL Lac objects (open circles) and FSRQs (closed circles) 

371 are quite distinct. The flux sensitivity (calculated in the same way as for the map shown in Fig. [2] 

372 and depicted as solid lines for two different galactic latitudes) is fairly strongly dependent on the 

373 photon index. The upper envelope in the spectral index - flux (>100 MeV) plot reflects that the 

374 peak sensitivity of the LAT is at energies much higher than 100 MeV. These ranges of spectral 

375 index and apparent flux limits translate to approximately constant limits above 1 GeV. For a 

376 photon index of 2.2, the 10 a flux sensitivity Fioo — 5x 10~^ ph cm~^ s~^, about 3 times lower 

377 than that of the Third EGRET catalog. 



'''http://fermi.gsfc.nasa.gov/ssc/data/policy/LAT_Monitored_Sources.html 
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378 4.2. Flux 

379 It makes sense to compare the LBAS fluxes wit h those reported in the Third EGRE T Catalog 



3B0 for the EGRET sample. As several analyses (e.g. iMiicke &: Pohll 120001 : iDermerl 120071 ) used the 
3B1 peak flux (maximum flux in all EGRET viewing periods), instead of the mean flux because of 

382 the fairly non-uniform coverage in the EGRET Catalog, comparisons will be performed both for 

383 the mean and peak flux distributions. For EGRET, both distributions are biased as observations 

384 were preferentially made of sources known to be highly variable in the gamma-ray band, and some 

385 of the observations were triggered by ToO requests when an object was brightly flaring in other 

386 wavebands. No such bias exists for the LAT. 

387 Fig. [8^ compares the mean flux distribution measured in the LBAS with that measured in the 

388 EGRET sample. The high-flux ends of these distributions look similar. This observation points 

389 to a nearly constant global gamma-ray luminosity of detectable blazars at a given time, as can 

390 naively be expected. In stark contrast, the weekly peak flux distributions (Fig. [8b) look different, 

391 the peak fluxes being significantly higher in the EGRET sample. This feature probably arises from 

392 the shorter sampling period for the Fermi-hAT as compared to EGRET. In the 3-month period 

393 considered here, a given source had much less opportunity to explore very different states than in 

394 the 4.5 years over which the EGRET observations were conducted. Another illustration of this 

395 effect is given in Fig. [Htjd where the peak flux vs the mean flux and the peak flux/mean flux ratio 

396 distributions are shown respectively. The inference that the gamma-ray blazars have characteristic 

397 variability timescales of months to years is well confirmed by the observation that only ~ 30% of 

398 the EGRET blazars are still detected by the LAT at a comparable flux. 



399 4.3. Photon index 

400 The photon index distribution gives insight into the emission and acceleration processes acting 

401 within the AGN jets, as it enables some of the physical parameters involved in these processes to 

402 be constrained. Moreover, it can be used to test whether the BL Lac and FSRQ populations have 

403 different 7— ray emission properties. 

404 Fig. [9] top displays the photon index distribution for a ll the LBAS sources. Th is distribution 



405 looks fairly similar to that observed for the EGRET sample (jNandikotkur et al.ll2007l ): it is roughly 

406 symmetric and centered at 7= 2.25. The corresponding distributions for FSRQs and BL Lacs are 

407 shown in Fig. [9] middle and bottom respectively. These distributions appear clearly distinct, with 

408 little overlap between them. This is a remarkable feature, given that the statistical uncertainty 

409 typically amounts to 0.1 for most sources. The distributions have (mean, rms)=(1.99, 0.22) for 

410 BL Lacs and (2.40, 0.17) for FSRQs. We used a Kolmogorov-Smirnov (KS) test to test the null 

411 hypothesis that both index samples are drawn from the same underlying distribution and found 
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412 a probability of 2 x 10~^^ Although indications for the existence of two spectrally distinct 

413 populations (BL Lacs and FSRQs) in the E GRET blazar sample were mentioned in the literature 

414 (|Pohl et alj 119971 : IVenters &: Pavlidoul 120071 ) . this is the first time that the distinction appears so 

415 clearly. The mean photon index of the 10 HBLs included in the LEAS is 1.76, i.e significantly lower 

416 (sources are harder) than the mean of the whole BL Lac subset as expected for these high-energy 

417 peaked sources. 

41B To infer physical properties of the blazar populations from the observed photon index distri- 

419 butions, possible instrumental and/or statistical effects have to be assessed. A systematic bias may 

420 indeed arise in the likelihood analysis of sources with low photon statistics. To quantify this possi- 

421 ble bias we performed a simulation study with the gtobssim tool which is part of the ScienceTools. 

422 This tool allows observations to be simulated using the instrument response functions and the real 

423 orbit/attitude parameters. Both instrumental and diffuse backgrounds were modeled on the basis 

424 of the real backgrounds observed by the LAT. 



425 1. Samples of sources (100 FSRQs and 100 BL Lacs) with random positions in the |6[ >10° sky 

426 were simulated. 

427 2. The real spacecraft orbit and attitude profiles spanning 94 days starting from Aug 4 2008 

428 were used. 

429 3. The sources were assumed to have a power-law energy distributions. The photon index was 

430 drawn from a gaussian distribution with (mean, sigma)= (2.0,0.3) for BL Lacs and (2.3,0.3) 

431 for FSRQs. These distributions are referred to as "input" probability distribution functions 

432 (pdfs). 

433 4. Fluxes were generated according to a lognormal distribution f{x) = exp with 

434 /i = In 10^'^ and a = 0.4 

435 5. A likelihood analysis was performed for all sources. The pdfs of the spectral indices and fluxes 

436 were built for sources with TS>100 ("like" pdfs). The TS cut was also applied to the "input" 

437 pdfs. 



438 Possible bias arising from the likelihood analysis as well as the robustness of the separation 

439 between BL Lac and FSRQ "like " pdfs were studied by means of KS tests. "Input" and "like" pdfs 

440 were found to be consistent with a probability of 99.5%, 88.4% for BLLacs and FSRQs respectively, 

441 excluding any sizeable bias coming from the likelihood analysis. The TS cut was observed to only 

442 affect the distribution tails. Concerning the separation between BL Lacs and FSRQs, the KS test 

443 returned that the probability for the two distributions to result from the same parent distribution 

444 is 7 X 10~^. 



^We are aware of the fact that the KS test is not optimal for binned data, but it is accurate enough to reject the 
null hypothesis 



445 



5. Sources already detected by EGRET 



446 After an elapsed time of about 10 years, it is interesting to look at the fraction of the AGNs 

447 that were active in the EGRET era and are detected again by the LAT with a comparable flux. 

448 Out of 116 sources in the Fermi-L AT sample, 3 sources have positions compatible with sources in 

449 the Third EGRET Catalog. Two additional sources, OFGL J1802.6-3939 and OFGL J1833.4-2106 

450 located at |b| <10° fulfills this condition as well. The 35 sources are listed in Table 4, along with 

451 the mean fluxes and photon indices measured by the Fermi-LAT and EGRET as well as the AGN 

452 class. These 35 AGNs are composed of 20 FSRQs, 11 BL Lacs, 3 of uncertain type and 1 AGN 

453 (Cen A). The BL Lacs are again overrepresented (with a fraction of 31%) as compared to the 

454 1st year sky survey EGRET sample (14 out of 60, i.e. 23%). The (non-simultaneous) fluxes and 

455 indices measured by both intruments are compared in Fig. [T2j The large scatter observed when 

456 comparing the fluxes (Fig. [12] left) can be expected from the variable nature of the blazar emission. 

457 The scatter observed when comparing the photon indices is more moderate, as could be expected 

458 from the fairly strong correlation between photon index and blazar class mentioned above. For 

459 many sources, and most especially for BL Lacs, the indices are measured by the Fermi-hAT with 

460 a much better accuracy. 



461 6. Radio gamma-ray connection 

462 With 116/125 high \b\, non-pulsar LAT bright sources associated with radio sources in the 

463 CRATES/CGRaBS and the Roma-JiZCAT lists, we confirm the findings of the 3EG catalog. In 

464 particular, 98/106 (~ 92%) of our high confidence associations have flux density above 100 mJy at 

465 8.4 GHz. In terms of the radio luminosity L,. = vL{v), the sources in the present sample with a 

466 measured redshift span the range lO^^-*^^ < Lf < 10^^'^^ ergs~^. As shown by the histogram in Fig. 
467 1131 BL Lacs and FSRQ are not uniformly distributed in this interval, with the former on average 

468 at lower radio luminosities (Log L^, bl Lacs = 42.8 ± 1.1 [ergs""*^]) than the latter (LogL^^FSRQ = 

469 44.4 lb 0.6 [ergs~^]). Blazars of uncertain type generally lack a redshift. Of the two radio galaxies 

470 associated with objects in the LBAS, NGC 1275 is similar to BL Lacs {L^ = 10^^'^^ ergs~^), while 

471 Cen A lies at the very lower end of the radio power distribution, with = 10^^'^^ ergs~^. 

472 Cen A, the source associated with OFGL J1325.4— 4303, is also the only source showing a 

473 significant amount of extended radio emission at low frequency (5'8.4/5'iow = 0.005). For all other 

474 sources with a low frequency (typically, 365 MHz from the Texas survey, 325 MHz from the WENSS, 

475 or 408 MHz from the B2) and a high frequency, high resolution (typically at 8.4 GHz from CRATES) 

476 fiux density measurement, we find little or no evidence of significant deviation from Liow = -^^8.4- 

477 Therefore, we find not only that all the sources in our sample are radio emitters, but that they 

478 also possess compact cores with flat radi o spectral index and mu ch higher luminosity than those of 



479 radio galaxies of similar or larger power (jGiovannini et al.lll988l ). 



480 Thanks to the comparatively large number of LBAS sources, it is worthwhile to perform a 
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481 statistical comparison of their properties in the gamma-ray and radio ba nds. Previous studi es based 



482 on EGRET data for 38 extragalactic point sources have been reported (jMiicke et al.lll997l ). which 

483 did not support claims of correlations between radio and gamma-ray luminosities. In particular, 

484 the analysis of possible correlations needs to be treated with care, because of the many biases that 

485 can arise, e.g. from the common redshift dependence when one considers luminosities, or from the 

486 reduced dynamical range when one considers mean flux densities, just to name a few. 

437 We have therefore looked at several possible pairs of observables, and we summarize our results 

438 in Table [5j In general, we apply the K-correction to the luminosities but not to the fluxes, since 

439 this would introduce a bias for the sources without a known redshift. We show in Fig. [T3] (left 

490 panel) the peak gamma-ray flux S'£;>iooMoV vs the radio flux density SsagHz from CRATES (or 

491 NED, in the few cases in which the source is not in the CRATES list). In general, BL Lacs tend 

492 to populate the low flux region, and FSRQs the high flux region. Such a constellation is prone 

493 to create correlations artificially from purely combining both populations. Given their different 

494 redshift distributions, this would be even more apparent in the luminosity plane. For this reason, 

495 it is necessary to consider the two populations separately (see Table [5]). Indeed, the results of our 

496 analysis show the significance of a radio-to-gamma-ray connection to be marginal at most on the 

497 basis of the present data, in particular for the FSRQs. Clearly, there is need for a deeper analysis 

498 on an enlarged sample regarding this issue, including Monte-Carlo simulations, which we defer to 

499 a forthcoming paper. 

500 Finally, we show in the right panel of Fig. [TH the radio luminosity vs. gamma-ray spectral 

501 index plane. Thanks to the large LAT energy range, the separation between BL Lacs and FSRQs is 

502 readily seen, showing a trend of softer spectral indices for more luminous radio sources. Moreover, 

503 this plot seems quite effective at finding sources of a different nature, such as the radio galaxy 

504 Cen A, whose gamma-ray index is much softer than t hat of other low p ower radio sources. For 



505 instance, OFGL J00174-0503 is a FSRQ at z = 0.227 (IHealev et al.ll2008l ) with index = 2.71 and 

506 radio luminosity Lr = 10^^'^^ ergs~^, which could then be a rare case of low-energy peak and low 

507 radio luminosity blazar. The other source with large photon index (2.60) and comparatively low 

508 radio luminosity {Lr = 10^^'^^ ergs"^) is associated with the peculiar source PMN J0948+0022. 



509 7. Population Studies 

510 As described before, the LBAS includes 57 FSRQs, 42 BL Lac objects, 5 blazars of uncertain 

511 type and 2 radio galaxies. Ten other sources have lower confidence associations with known blazars. 

512 This sample is already comparable with that provided by EGRET and can be used to derive some 

513 early results about the redshift and source count distributions and the luminosity function of blazars. 
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514 7.1. Redshifts 

515 Fig. [15] and Fig. [16] display the redshift distributions for FSRQs and BL Lac objects, respec- 

516 tively, and their comparison with those of the parent distributions in the BZCat catalog. Please 

517 note that 12 of 42 BL Lacs have no measured redshifts. BL Lac objects are generally found at low, 

518 z < 0.5, redshift, whereas the peak of the FSRQ r edshift distribution is a round z = 1. Similar dis- 



519 tributions were observed for the £^Gi?i?T blazars (iMukherjee et al .11 199 71 ). In the future, as fainter 

520 sources become visible, detection of additional nearby radio-galaxies will enhance the number of 

521 very low redshift objects in the AGN redshift distributions measured with the Fermi LAT. 

Fig. [13 shows the luminosities of the detected sources plotted as a function of their redshifts. 
The isotropic gamma-ray luminosity was derived using: 

= 4^54/(1 + (1) 

522 Here, S is the 7-ray energy flux (E > 100 MeV), a is the energy index and d^, is the luminosity 

523 distance. A beaming factor 5 = 1 was assumed. The solid curve corresponds to a flux limit of Fioo 

524 = 4 X 10~* ph cm~^s~-'^. 



525 7.2. log N - log S 

526 7.2.1. Monte Carlo Simulations 

527 Proper population studies must rely on a thorough understanding of the properties of the 

528 survey where these objects have been detected. In order to properly estimate the source-detection 

529 efficiency and biases, we performed detail ed Monte Carlo simu lations. The me thod we adopted 

530 is the one developed for ROSAT analysis (jl lasinger et al. 1993 ) and lately used ( Cappelluti et al 



531 120071 ) for the analysis of the XMM-COSMOS data. For each source population (blazars, FSRQs 

532 and BL Lacs) we created a set of >20 LAT all-sky images with background patterns resembling 

533 as close as possible the observed ones. The simulations were performed using a similar method as 

534 that described in section 4.3. An extragalactic population of pointlike sources was added to each 

535 simulated observation. The coordinates of each source were randomly drawn in order to produce 

536 an isotropic distribution on the sky. Source fluxes were randomly drawn from a standard log A^- 

537 log S distribution with parameters similar to the one observed by LAT (see next section) . Even 

538 though the method we adopt to derive the survey sensitivity does not depend on the normalization 

539 or the slope of the input log A^-log S", using the real distribution allows us to produce simulated 

540 observations which closely resemble the real LAT sky. The photon index of each source was also 

541 drawn from a Gaussian distribution with observed mean and 1 a width consistent with the real 

542 population. Thus, for the three simulation sets we adopted the following photon indices similar to 

543 the measured ones: 



544 • 2.24ib0.25 for the total blazar population; 
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545 • 2.41±0.17 for the FSRQ population; 

546 • 1.98ib0.22 for the BL Lac population. 

More than 30000 sources were simulated for each population. The mock observations were 
processed applying the same filtering criteria used for real in-flight data. Source detection was 
performed on E>200MeV photons with a simplified version of the detection algorithn|l|. For 
every pair of input-output sources, we computed the quantity: 

547 where x, y and 5* are the source coordinates and flux of the detected sources while xq, yo and Sq are 
54B the corresponding values of the input sources and ax, (Ty, as the associated statistical uncertainties. 

549 We then flagged those with the minimum value of as the most likely associations. Only sources 

550 at lb[> 10° are retained. 

551 The goal of these simulations is to derive the probability of detecting a source (with given 

552 mean properties, e.g. photon index and flux) in the LAT survey as a function of source flux. This 

553 can be computed from the simulations reported above as the ratio between the number of detected 

554 and input sources in a given flux bin. The detection efficiencies for the three source populations 

555 are reported in Fig. [TSl A few things can be noted readily. First, the bias of the LAT survey 

556 against soft sources (i.e. FSRQs) is apparent. Second, the LAT |b| >10° survey becomes complete 

557 for F(>100 MeV)>2x 10~^ ph cm~^ s~^ irrespective of the source photon index or its location in 

558 the sky. Multiplying these functions by the solid angle J7 of the survey (34089.45 deg^ in case of a 

559 |b|> 10° cut) yields the so called sky coverage which is used for the statistical studies reported in 

560 the next sections. 



561 7.2.2. Incompleteness of the Extragalactic Sample 

562 We report in Table [6] the composition of the |b| >10° sample. The number of sources with 

563 high-confidence associations is 106. Of these 57 are FSRQs and 42 are BL Lacs. As already shown 

564 in the previous sections, FSRQs and BL Lacs are represented in almost equal fractions in the LAT 

565 survey. The 5 blazars with uncertain classifications are likely split between these two categories 

566 as the redshift-luminosity plane (Fig. [T71) shows. The incompleteness factor varies as a function of 

567 the sample under study. When considering the non-pulsar part of the high-confidence sample, the 

568 incompleteness is given by low confidence and unassociated objects. This turns out to be ~11 %. 

569 However, when considering the FSRQ and BL Lac samples separately one must also include the 



^The complexity of the official detection algorithm makes it virtually impossible to apply it to a large number of 
data sets. We tested on real data that, for the scope of this investigation, our simplified detection algorithm produces 
results consistent with more elaborate ones. 
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570 sources with uncertain classifications. Thus the incompleteness factor of the FSRQ and BL Lac 

571 samples rises to ~15 %. A reasonable and simple hypothesis is one which assumes that these sources 

572 reflect the composition of the identified portion of the sample. This would mean that there are 

573 an additional ~9 FSRQs and ~7 BL Lacs are hiding among the unidentified/unassociated/low- 

574 confidence sources. These uncertainties will be used in the next sections. 

575 Since the uncertainty due to the incompleteness is relatively large, we will use a flux-limited 

576 sample to verify the results derived from the main sample. Indeed, for Fioo > 1.25 x 10^^ ph cm~^ 

577 s~^ the number of uncertain, unassociated, and low-confidence sources falls to 2 (and 2 are anti- 

578 associated). Above this flux limit, the sample contains 44 sources of which 29 and 9 are FSRQs and 

579 BL Lacs respectively, while 2 are Radio galaxies. Moreover, all but one BL Lac have a measured 

580 redshift. Thus, while low numbers penalize this flux-limited sample, its incompleteness is <5%. 



581 7.2.3. Source Counts Distributions 

The source counts distribution, also known as the log A'^-log S, flux, or size distribution, is 
readily computed once the sky coverage is known through the expression: 

Ni>S) = J2^,^^S~', (3) 

i=l * 

where Ng is the total number of detected sources with fluxes greater than 5*, and fij (i.e. Fig. \TE\ 
multiplied by the solid angle) is solid angle associated with the flux of the ith source. The variance 
of the source number counts is defined as 

Ns . . .2 

582 In building the source counts distributions, we used the source flux averaged over the three month 

583 timescale. The log A^-log S of the entire extragalactic sample (excluding pulsars) is shown in 

584 Fig. m 



We fitted the source counts distribution with a power-law model of the type: 



A common practice in this case (e.g., see lAiello et al.ll2008l ) is to fit the unbinned dataset employing 



a maximum likelihood (ML) algorithm. For this purpose the ML estimator can be written as 

r^-2Eln "'fg<^-' , (6) 

i J n{S)n{S)dS 

585 where i runs over the detected sources. The 1 a error associated to the fitted parameters (in this 

586 case a) is computed by varying the parameter of interest, while the others are allowed to float. 
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until an increment of A£ =l is achiev ed. This gives an estimate of the 68 % confidence region for 



588 the parameter of interest (|Avnilll976l ). In this formulation of the ML function, the normalization 

589 A is not a parameter of the problem. Once the slope a is determined, the normalization is derived 

590 as the value which reproduces the number of observed sources. An estimate of its statistical error 

591 is given by the Poisson error of sources used to build the log A'"-log S. 

592 Since the sky coverage is somewhat uncertain at very low fluxes, the fit is performed above 

593 F(>100 MeV)= 7 x 10~® ph cm~^ s~^ even though all the data are displayed. The results of the 

594 best fits to the different source counts distributions are summarized in Table [71 It is clear that all 

595 distributions are compatible, within their errors, with a Euclidean distribution (a = 2.5). In order 

596 to check the stability of our results we have shifted the sky coverage of Fig. [18] by 20 % on either 

597 side. Taking the whole extragalactic population as an example (see first line of Table [7]) we get that 

598 the best fit values of the slope are 2.47 and 2.62 for the —20 % and +20 % case respectively. These 

599 values are consistent within the error (e.g. 2.59+0.12), showing that at bright fluxes our analysis 

600 does not suffer from major systematic uncertainties in the sky coverage. The same result holds for 

601 the other log A^-log S distributions reported in Table \7\ 

602 The log A^- logS distributions for FSRQs and BL Lacs are shown in Fig. [20] and [2TJ We do 

603 not find any indication of a break in the source counts distributions of the two populations. As 

604 the fitting results of Table [7| show, there might be an intrinsic difference between the log A^-log S 

605 of both populations, with the source counts distribution of BL Lacs being flatter than that of of 

606 FSRQs. However, both of them are compatible within 1 a errors with the Euclidean value of 2.5. 

607 Moreover, the analysis of the flux-limited sample (see bottom part of Table [7]) conflrms the results 

608 of the main sample, showing that incompleteness is not a main issue in this study. 

609 For the EGRET sample, a surface density for Finn ^ 10~''ph cm~^ s^^ of FSRQs and BL Lacs 



610 of 3.31 sr^^ and 0.83 sr^^, respectively, is reported ( Miicke Sz Pohlll200d ). From LAT we derive that 

611 the surface density (above Fmn > 10"^ ph cm~^ s~^) of FSRQs and BL Lacs is 4.41+0.72 sr~^ and 

612 1.01+0.17 sr~^ respectively. Thus the LAT results are in good agreement with EGRET. 

A measurement of the number fluence using the average three-month fluxes of bright Fermi 
blazars of different classes is readily obtained from the log A^-log S distributions through the 
expression: 

rfmax 

^diffuse = I -j^ S dS . (7) 



rnin 



613 Unless otherwise stated, we adopt a value for fmin of 4x10 ph cm ^ s ^. To compare with the 

614 EGRET results, the upper limit of integration cannot be set to inflnity. Indeed, all point sources 



615 detec ted above Fmn ~ 10 ''ph cm ^ s ^ in the Second EGRET Catalog (2EG; [Thompson et al 



6i6 ll995ll were subtracted in the measurement of the extragalactic diffuse 7-ray background (EDGB) 
617 (jSreekumar et al.lll998l ). Thus, we set fmax to 10~^ph cm~^ s~^. The integral in Eq. [7] yields a 



618 total flux of 1.06(+0.09) X 10 ^ ph cm ^ s ^ sr ^ This can be compared with the intensity of the 

619 EDGB, as measured by EGRET, of 1.45x10"^ ph cm~^ s~^ sr~^. Already in this small flux range. 
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620 LAT is resolving into pointlike sources ~7% of the EGRET FjT)GH. Preliminary analysis of the log 

621 A^-log S distributions shows that LAT is expected to resolve a much larger fraction of the EDGB 

622 within the next few months of observation. 



7.3. Evolution of Blazars 

7.3.1. Evolutionary Test 



625 A simple and robust test of evolution is the V/Vmax test ISchmidtl (|l968l ). The quantity 

626 V/Vmax is the ratio between the (comoving) volume within which the source has been detected 

627 and the maximum comoving volume available for its detection. For a given source, V/Vmax is 

628 expected to be uniformly distributed between and 1. For a population uniformly distributed 

629 in Euclidean space (and with constant properties with z) and non-evolving, the average V/Vmax 

630 should be consistent with a value of 0.5. The error on the average value is c = l/(12A^)i/2 for N 

631 sources. A value of <V/Vmax> > 0.5 indicates positive evolution (more sources or brighter sources 

632 at earlier times), and the opposite indicates negative evolution. 

The comoving volume for a ith source is given by 

l-z=z, 

V= — n{L„z)dz, (8) 

./^=n clz 



633 wher e dV/dz is the comoving volume element per unit redshift and unit solid angle (see e.g. iHogg 



19991 ) and r2(Lj,z) is the aforementioned sky coverage for the source with rest-frame luminosity 



635 Li at redshift z. We note that t he definition of the W Vmax reported in Eq. [8] encompasses also 



636 the definition of the Ve/Va test (jAvni &: Bahcalll Il980l ). which for the purposes here are formally 

637 equivalent. 

638 We computed the average <V /Vmax> for FSRQs, BL Lacs and all sources in the high- 

639 confidence sample with measured redshift (these includes the sources with uncertain classification). 

640 The results are summarized in Table [H All 57 FSRQs present in the extragalactic sample (see 

641 Table [6|) have a measured redshift. The V/Vmax shows that the population of FSRQs detected 

642 by LAT evolves positively (i.e. there were more FSRQs in the past or they were more luminous) 

643 at the 3 a level. This result is also confirmed by the analysis of the 29 FSRQs which constitute a 

644 flux-limited sample (see lower part of Table [6|). 

645 Only 31 out of the 42 BL Lac objects have a measured redshift. The V/Vmax test is compatible 

646 within ~ la with no evolution. Assigning the mean redshift value of the BL Lac sample (i.e. 

647 (z)=0.38) to those objects without a without a redshift produces a value of (y/yMAx)=0.472ib0.046. 

648 The result does not change if the redshift is drawn from a Gaussian distribution with mean and 

649 dispersion consistent with the observed redshift distribution of BL Lacs. However, it is difficult to 

650 assess the validity of both these hypotheses. Indeed, the fact that these objects sh ow a featureless 



651 continuum might suggest that their redshifts could be the largest in the sample (jPadovani et al 
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20071 ). In this case, their true redshift would produce a larger value of the V/Vmax statistic. The 



653 V/Vmax of all the objects with a measured redshift in the high-confidence sample is compatible 

654 with no evolution. 



7.3.2. Luminosity Function of FSRQs 



We estimate the gamma-ray luminosity function (GLF) in fixed redshift bins using the 1/Vmax 
method (equivalent in our formalism to the 1/Va method). For each bin of redshift the GLF can 
be expressed as 

dN 1 ^ 1 



A 



Vmax^i 



656 where VMAX,i is the maximum comoving volume associated with the ith source (see Eq. [8|) . The 

657 cumulative and differential luminosity functions of FSRQs, in three redshift bins, are reported in 

658 Fig. [22l One thing is readily apparent from this figure. FSRQs are strongly evolving. A non- 
659 evolving population would have GLFs which are continuous across different redshift bins. In the 

660 case of FSRQs we note a change in space density (or luminosity) with redshift. Also, one can see 

661 that the space density of intermediate-luminosity FSRQs (e.g. Ly ~10'^''erg s~^) is increasing with 

662 redshift. On the other hand, the most luminous FSRQs have an almost constant space density 

663 with redshift. This might be a sign of a cut-off in the evolution of FSRQs. A decline ii i the space 



664 density of luminous FSRQ s has also been det ermined at radio and X-ray energies (e.g., IWall et al 



665 12OO5I : iPadovani et al.l 120071 : lAiello et al.l |2009| ) . We derive from the GLF that the space density of 
666 FSRQs with L^ > 7 x 10"^^ erg s'^ is 1.05±0.13 Gpc'^ 

We made a Maximum Likelihood fit to the three unbinned datasets using a simple GLF model 
defined as 



$(L 



7; • 



The ML estimator can be expressed similarly to Eq. [6] by the expression 



-2Vln 



1 Zi)V {Ly^i, Zi) 



(10) 



(11) 



667 The results of the ML fits to the GLF of FSRQs are summarized in Table [H For z > 1 , the GLF 

668 can be successfully parametrized by a single power-law model. The slope is compati ble with the 



669 canonical value of 2.5-2.8 determined for X-ray selected samples of radio-quiet AGNs (jUeda et al 



2003 



Hasinger et al 



2005 



Silverman et al.ll2008l ). This indicates that at high redshifts the Fermi- 



671 LAT is sampling the bright end of the luminosity distribution of FSRQs. For z < 1, the best-fit 

672 value of the slope /3 is 1.56=fc0. 10, comp atible with luminosity function slopes found in radio/X-ray 

673 selected samples (jPadovani et al.l 120071 ). This is much flatter than the canonical value of /3 = 2.5. 

674 As the cumulative GLF shows (left panel of Fig. [22]), there might be a hint of a break with respect 

675 to a simple power-law model in the GLF. A more detailed analysis, comparing different methods 

676 to derive the GLF and its evolution, will be considered in future publications. 
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7. 3. 3. Luminosity Function of BL Lacs 



678 The luminosity function of BL Lacs, reported in Fig. [231 is in agreement with the results of 

679 the V/Vmax test. Indeed, sub-dividing the entire BL Lac sample in two bins of redshift produces 

680 two GLFs which connect smoothly to each other. A simple power-law GLF describes the entire 

681 dataset well. The GLF slope is (3 = 2.17 ± 0.05 and is well in agreement with t he value of 



682 2.12zb0.16 reported for a radio/X-ray selected sar nple of BL Lacs (IPac 



683 of 12 EGRET BL Lac objects in a recent study (jBhattacharva et al. 



ovani et al.ll2007l 'l. The GLF 



20081) was found to show no 



684 signi f icant evidence for evolu tion, with a GLF slope /3 = 2.37 ±0.3. Past claims (e.g.. iRector et al 



2000; iBeckmann et al.ll2003l ) of negative evolution of BL Lac objects, selected mainly in the X-ray 



686 band, are not confirmed by our data. The dynamical range of the LAT GLF samples 4 decades 

687 in luminosity and nearly 8 in space density. From our GLF we derive that the density of BL Lac 

688 objects with L^ >3xl0'^^erg is 1.9(±0.4) x 10""^ Mpc^^. 

689 Above a luminosity of L^ ~ lO'^'^erg s~^, the cumulative density of BL Lacs and FSRQs is 

690 comparable, with BL Lacs being ^ 3 times less numerous than FSRQs. However, given the fact 

691 that they reach lower luminosities, BL Lacs are ^ 200 times more abundant than FSRQs above 

692 their respective limiting luminosities. 



8. Discussion 



694 The value TS > 100 defining the detection significance for bright sources corresponds to > 10 u 

695 significance, or a limiting flux over the entire high-latitude sky of ~ (3 - 10) x 10~® ph(> 100 MeV) 

696 cm^^ s^^ during the three-month sky survey. In comparison, EGRET reached a 5a high-confidence 

697 on-axis flux limit of ~ 15 x 10^*^ ph(> 100 MeV) cm~^ s~^ for a two-week pointing over ~ 0.5 sr 



698 of the sky, only becoming complete at S -Fioo~ 25 x 10 ^ ph (cm ^ s ^ ( Dermer 2007). Of the 



699 66 hi gh-confidence and 27 lower-confidence AGN associations in the 3EG catalog (jHartman et al 



19991 ). 32 sources in the Fermi-LAT sample were also detected with EGRET. An additional source 



701 is detected at |b|< 10°. Many of the other high-confidence EGRET sources are detected with 

702 Fermi-LAT at TS < 100, reflecting the rapid variability and periods of activity of 7-ray blazars on 

703 timescales of years or longer. 

704 During the 18-month EGRET all-sky survey when exposure to all parts of the sky was relatively 

705 uniform compared to the remain der of the mission, 6 high-confidence blazars consisting of 14 BL 

706 Lacs and 46 FSRQs were found (jFichtel et al.lll994l ). Compared with ~ 23% of EGRET blazars 

707 being BL Lac objects, nearly 40% of the Fermi-LAT blazars are BL Lac objects. The larger fraction 

708 of BL Lac objects in the Fermi bright AGN sample is partly a consequence of the good sensitivity 

709 to high-energy emission by Fer mi-LAT , whereas self- ve toing in EGRET reduced its effective area to 

710 photons with energies > 5 GeV (jThompson et al.lll993l ). Consequently, dim hard-spectrum sources 

711 are favored to be detected with the Fermi-LAT compared to EGRET. 
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712 A clear separation between the spectral indices of FSRQs and BL Lacs is found in the Fermi- 

713 LAT data (Fig. [9]), with mean photon indices of T = 2.40±0.17 (rms) for FSRQs and F = 1.99±0.22 

714 (rms) for BL Lac objects. A KS test gives a probability of 2 x 10~^^ for the two index samples to be 

715 drawn from the same parent distribution. Moreover, the SEDs of bright flaring blazars in the cases 

716 of 3C 454.3 and AO 0235+164 show a spectral softening at E > 2 GeV. If this behavior persists 

717 in weaker FSRQs, then an even greater fraction of BL Lac objects will be found in Fermi-hAT 

718 analyses over longer times, because signal-to-noise detection significance for weak hard-spectrum 

719 sources becomes better than for weak soft-spectrum sources due to the reduced background at 

720 higher photon energies. 

721 Another reason for the larger fraction of BL Lac objects in the Fermi-hAT blazars could be 

722 related to the redshift distribution of the bright AGNs. The BL Lac objects are dominated by 

723 low-redshift, z < 0.5 blazars, with a tail extending to z ~ 1, whereas the FSRQs have a broad 

724 distribution peaking at z Rr: 1 and extending to z ^ 3 (see Fig.fTSland Fig . [THll . These distributions 



725 are similar to the distribution of EGRET hlazais (iMukherjee et al.lll997l ). Because the peak of the 

726 EGRET FSRQ redshift distribution is already at z « 1, detection of higher redshift FSRQs with 

727 the more sensitive Fermz-LAT would be impeded by cosmological factors that strongly reduce the 

728 received fluxes. Moreover, the period of dominant AGN activity was probably at z ~ 1 or 2. The 

729 increased sensitivity for the BL Lac objects with Fermi-LAT, on the other hand, allows it to probe 

730 beyond the low-redshift population of BL Lac objects detected with EGRET where the detectable 

731 volume is still rapidly increasing with z. The likelihood of detecting z ~ 1 BL Lac objects does 

732 depend, however, on their evolution. 

733 The simplest index of population evolution is the V/Vmax test. We found (V/Vmax) = 

734 0.43 lb 0.055 for the BL Lac objects with redshift in the LBAS (Table 8), so that the BL Lac objects 

735 are within ~ la of showing no evidence for evolution. For the FSRQs in the LBAS, by contrast, we 

736 found (V/Vmax) = 0.64 it 0.04, so that the FSRQs exhibit strong positive evolution. The strong 

737 positive evolution of FSRQs and weakly negative or no evolution of BL Lac objects in the LBAS is 

738 contrary, however, to our reasoning that population evolution of the lower redshift BL Lac objects 

739 explains the larger fraction of BL Lacs in the LBAS compared with the BL Lac fraction observed 

740 with EGRET. As indicated by the indices of the log - log 5" (eq. 6 and Table 7), which show 

741 much weaker evidence for evolution than given by the V/Vmax test, the actual situation may be 

742 more complicated and depend on both density and luminosity evolution. 

743 The BL Lac objects are found to display systematically harder spectra, with vF^, spectra rising 

744 at GeV energies, compared to the powerful FSRQs where the peak of the uF^, spectrum is at photon 

745 energies < 100 MeV - GeV. This is generally attributed to a different dominant radiation process; 

746 self-Compton scattering of the jet electron's synchrotron emission in the case of BL Lac objects, 

747 and Compton scattering of external radiation field s in the case o f FSRQs if leptonic processes 

748 dominate the radiation output (recently reviewed in iBottcherlboO?! ). The excellent sensitivity and 

749 full-sky coverage of the Fermi LAT is, for the first time, giving us broadband evolving SEDs from 

750 the radio to the 7-ray regime in sources like 3C 454.3, PKS 2155-304, and others that will require 
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751 detailed spectral modeling to assess the relative importance of self-Compton and external Compton 

752 scattering processes in the different blazar classes. 

753 Such results will be important to determine whether FSRQs and BL Lac objects may have 

754 a direct evolutionary relationship, or instead represent separate unrelated tracks of supermassive 

755 black hole fueling and growth. A scenario whereby BL Lac objects are the late stages of FSRQs, 

756 as th e gas and dust produced i n a galaxy merger or tida l interaction fuels the supermassive black 

757 hole (jBottcher &: Dermeiil2002l : ICavaliere &: D'Eliall2002l ). provides a framework to understand the 

758 blazar phenomenology and makes definite predictions about the relative black hole masses in the 

759 two classes. The more abundant scattered radiation and fuelin g in the evolution from FSRQ to BL 



760 Lac o bject would then lead to a blazar sequence like behavior (jFossati et al.lll998l : iGhisellini et al 



19981 ) if the amount of accreting matter controls black hole power and the surrounding radiation 



762 field. 

763 It is sti ll premature to compare the number of blazars in this bright source list with prelaunch 



764 predictions dM ucke fc Poh]||200i 



np 

3; 



Stecker fc Salamonll200ll : iNarumoto Totanill2006l : lDermeijl2007l : 
Inoue &: Totani 2008) made on the basis of differing assumptions, to sensitivities ~ 5a rather than 



766 1 0(7, and over different spans of time. Nevertheless, nearly complete surveys with far more sources 

767 than detected with EGRET are now available for calculating luminosity and number evolution, 

768 with implications that can be compared with results from the EGRET era. 

769 This study can be used to examine the observational basis for assuming an underlying radio/7- 

770 ray connection used to cal c ulate the blazar contributio n to the 7-ray background (jStecker &: Salamon 



19961 : iGiommi et al.ll2006l : INarumoto &: Totanil 120061 ). Figure O shows that except for a (at most) 



772 weak correlation of the brightest 7-ray blazars with the most radio-bright blazars, the 7-ray and 

773 radio fluxes display a large amount of scatter. Whether a stronger correlation can be found by 

774 comparing mean 7-ray fluxes with radio fluxes will require further study. But even at this early 

775 stage of the Fermi mission, we find that the bright sources can alre ady comprise about 7% of the 

776 diffuse extragalactic 7-ray background flux measured with EGRET (jSreekumar et al.lll998l ). 



777 We conclude this study by noting that the Fermi-L AT results imply the non-thermal luminosity 

778 density of AGNs on various size scales. A 7-ray blazar makes a contribution to the non-thermal 

779 emissivity oc L/V in terms of 7-ray luminosity and injection volume V derived from redshift. The 

780 Fermi-LAT results from Table 2 show that BL Lac objects provide local emissivities Ibl ^ 10^^ W 

781 Mpc~'^, whereas FSRQs have ipsRQ ~ lO^'^ W Mpc~^. Gen A, because of its proxir nity at d = 3.5 

782 Mpc , dominates the non-thermal luminosity, with £cenA ~ 3 x 10'^^ W Mpc~^ ( Dermer et al 



783 l2008l ). Sources of UHECRs must have a luminosity density within th e GZK radius, ~ 100 Mpc , 

784 of luHECR ~ 3 X 10^9 W Mpc~3 or Iuhecr ~ 10^'' ergs Mpc"^ yr"^ (|Waxman fc Bahcalll [l999l ) . 

785 To have sufficient emissivity within the GZK radius, if AGNs are the sources of the UHEGRs 

786 (jThe Pierre AUGER Collaboration et al.l 120071 ) , the Fermi-LAT results would therefore seem to 



787 favor BL Lac objects over FSRQs as the source of the UHECRs. 
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9. Summary 



7B9 We have presented a l ist of 116 bright, > 10a sources at |6| > 10° taken from the Ust of bright 

790 sources ()Abdo et al.ll2009d ) observed with the Fermi-LAT in its initial three-month observing period 

791 extending from August 4 to October 30 of 2008. Of these sources, 106 are associated with blazars 

792 with high confidence and compose the LBAS. The number of low-confidence AGN associations is 

793 11 (one source having two possible associations - one high and one low confidence). At |5| > 10°, 

794 5 sources out of a total of 125 non-pulsar sources remain unidentified. Two of the AGNs are 

795 associated with radio galaxies. The purpose of this work is to present the key properties of the 

796 AGN population of this bright GeV source list. The main results are summarized as follows: 



797 1. With a ~ 90% success rate from correlating the bright gamma-ray source list with AGN radio 

798 catalogs (CRATES/CGRABS, BZCAT) the bright extragalactic gamma-ray sky continues to 

799 be dominated by radio-bright AGNs. 

800 2. The number of HBLs in the LBAS detected at GeV energies (even when not flaring) has 

801 risen to at least 10 (out of 42 BL Lacs) as compared to one (out of 14 BL Lacs) detected by 

802 EGRET. Seven LBAS HBLs are known TeV-blazars. 

803 3. Only ~ 30% of the bright Fermi AGN list were also detected by EGRET. This may be a 

804 consequence of the duty cycle and variability behavior of GeV blazars. 

805 4. BL Lac objects make up almost half of the bright Fermi AGN sample (consisting of 57 FSRQs, 

806 42 BL Lac objects, 2 radio galaxies, and only 5 AGN remain unclassified), while the BL Lac 

807 fraction in the 3EG catalog was only ~ 23%. This feature most probably arises from the 

808 different instrument responses of the LAT and EGRET. 

809 5. The mean flux distribution of the Fermi AGN remains similar to the corresponding one based 

810 on the EGRET sample, while the peak flux distributions differ appreciably. 

811 6. We find a spectral separation between BL Lacs and FSRQs in the GeV gamma-ray band 

812 with FSRQs having significantly softer spectra than BL Lac objects. This confirms earlier 

813 indications for the existence of spectrally distinct populations in the EGRET blazar sample. 

814 The average photon index is 1.99ib0.22 (rms) for BL Lacs, with a tendency of HBLs displaying 

815 even harder spectra, and 2.40 ± 0.17 (rms) for FSRQs. A KS test gives a probability of 2 x 

816 10~^^ for the two index samples to be drawn from the same parent distribution. 

817 7. Fermi FSRQs in the bright source list are on average more luminous and more distant than the 

818 Fermi-detected BL Lac objects in that list. I.e., FSRQs exhibit a broad redshift distribution, 

819 starting with z = 0.158 (3C 273), peaking at z ~ 1 and extending up to z ~ 3 while BL Lacs 

820 are mostly found in the ~ 0.1 redshift bin with a tail extending up to z ~ 1. No significant 

821 relation between the gamma-ray photon index and redshift is found within each source class, 

822 in agreement with corresponding studies based on the EGRET AGN samples. 
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8. The peak gamma-ray flux is at best only weakly related to the 8.4 GHz radio flux, with the 
brightest gamma-ray AGNs having the largest radio flux densities. 

9. Using mean fluxes the Log N-Log S distribution of all the bright sources (except the pulsars) 
appears compatible with an Euclidean distribution without any breaks. This is also true 
within 1(7 for the source counts distributions of the FSRQ and BL Lac sample separately. 
Surface densities of 4.28 ± 0.72 sr^^ and 1.01 ± 0.17 sr'^ {Fwo > lO-'^ph cm'^ s"^) for 
FSRQs and BL Lacs, respectively, are reached. 

10. The combined emission in the flux range i^ioo.mcan ~ (7 - 10) X lO^^ph cm ^ s ^ observed 

831 from these individually resolved AGN during this three-month period already corresponds to 

832 ^ 7% of the EGRET detected extragalactic diffuse gamma-ray background. 

833 11. A V/Vmax analysis shows positive evolution at the 3 a level for the bright Fermi-detected 

834 FSRQs with the most luminous FSRQs having an almost constant space density with redshift, 

835 while for the Fermi-dctcctcd BL Lacs no evolution within one a is apparent. 

836 12. The gamma-ray luminosity function of bright FSRQs can be described by a single power-law 

837 with index ~ 2.5 and ~ 1.5 for the high (> 0.9) and low (< 0.9) redshift range, respectively, 

838 while the BL Lac gamma-ray luminosity function follows a power law with index ~ 2.1. The 

839 space density of gamma-ray emitting BL Lacs of ^ 190 Gpc~^ above their limiting luminosity, 

840 ~ 3 X 10^^ erg s~^, is a factor ^ 200 larger than for the Fermi-detected FSRQ population above 

841 their limiting luminosity, ~ 7 x 10^^ erg s~^. Thus, within the Fermi bright AGN list BL Lacs 

842 are intrinsically more numerous than FSRQs. Bright Fermi detected BL Lacs and FSRQs 

843 display comparable cumulative number counts above ~ lO^^erg s~^, with BL Lacs being ~ 3 

844 times more numerous than FSRQs. 

845 These early results from the first three months of the science mission of the Fermi Gamma ray 

846 Space Telescope demonstrate its exceptional capabilities to provide important new knowledge about 

847 7-ray emission from active galactic nuclei and blazars. As the Fermi-LAT data accumulate, many 

848 more AGNs at lower flux levels will likely be detected- as well as flaring AGNs at brighter fluxes 

849 than yet observed - helping to refine these results and improve our understanding of supermassive 

850 black holes. 
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Fig. 1. — Examples of weekly light curves for five bright blazars detected by Fermi-L AT and the 
Vela light curve for comparison (flux unit: xlO~^ photons cm~^ s— 1, please note the different 
scales). The dashed line is the average value and the grey area shows the 3% systematic error we 
have adopted. Different flux variability amplitudes and timescales are clearly visible in the blazar 
light curves. 
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Fig. 2. — Flux limit [E>100 MeV] (ph cm ^) as a function of sky location (in galactic coordi- 
nates), for a photon index=2.2 
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Fig. 3. — 95% error radius as a function of TS for the sources presented in this paper. FSRQs: 
closed circles, BL Lacs: open circles, Uncertain type: closed triangles, Radio galaxies: open stars. 




Fig. 4. — Location of the LBAS sources. FSRQs: closed circles, BL Lacs: open circles, Uncertain 
type: closed triangles, RG: open stars. 
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angular separation (sigma) 

Fig. 5. — Normalized angular separation between the Fermi-L AT location and that of the coun- 
terpart. The solid (dashed) histogram corresponds to the sources with high-(low-) confidence as- 
sociations. The solid curve corresponds to the expected distribution (x^ distribution with 2 d.o.f.) 
for real associations, the dashed one for accidental associations. 
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Fig. 6.— The X-ray (0.5 - 2.0 keV) vs radio flux density (1.4 GHz) plot of all X-ray detected 
blazars in the BZCAT catalog (small dots) and the Fermi-LAT detected blazars (BL Lacs: open 
circles, FSRQs: filled circles, blazars of uncertain type: triangles, radio galaxies: stars). 
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Fig. 7. — Flux [E>100 MeV] vs photon index for the 116 sources. FSRQs: closed circles, BL 
Lacs: open circles, Uncertain type: closed triangles, Radio galaxies: open stars. The solid curves 
represent the TS=100 limit estimated for two galactic latitudes b=20° and b=80° . The dashed 
curve represents the TS=100 limit for b=80° and 0.2 <E< 3 GeV. 
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Fig. 8. — a) Comparison of mean flux distribution for blazars detected by Fermi-LAT (solid) and 
EGRET (dashed), b) same of as a), for the peak flux distribution, c) Peak flux as a function of 
mean flux, for the Fermi-LAT (closed circles) and EGRET (open circles) AGNs. d) same as a), for 
the peak/mean flux ratio. 
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Fig. 9. — Photon index distributions for the LBAS blazars. Top: Ah sources. Middle: FSRQs. 
Bottom: BL Lacs. 
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Fig. 10. — Gamma-ray SED of 3 bright blazars calculated in five energy bands, compared with 
the power law fitted over the whole energy range. Left: 3C454.3 (FSRQ), middle: AO 0235+164 
(IBL), right: Mkn 501 (HBL) 
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Fig. 11. — Left: LBAS photon index as a function of redshift. Same symbols as before. Right: 
same as left, for the EGRET sample. 




Fig. 12. — Left: Fermi-LAT vs EGRET mean flux for the 33 AGNs present in both samples 
FSRQs: closed circles, BL Lacs: open circles, Uncertain type: closed triangles. Right: same as left, 
for photon index. 
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Fig. 13. — Histogram of the radio power distribution for LBAS sources, for all sources (upper 
panel), FSRQs (middle), and BL Lacs (bottom) only. 
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Fig. 14. — Radio vs. gamma-ray properties. Left: peak gamma-ray flux vs. radio flux density at 
8.4 GHz; the dashed hues show the CRATES flux density Umit and the typical LAT detection 
threshold. Right: gamma-ray photon index vs. radio luminosity. 
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Fig. 15. — Redshift distribution for the FSRQs in the LBAS (soHd) and in the BZCat catalog 
(dashed). 
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Fig. 16. — Redshift distribution for the BL Lacs in the LBAS (sohd) and in the BZCat catalog 
(dashed). 
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Fig. 17. — Gamma-ray Luminosity vs redshift for the LB AS. The soUd line was drawn using a 
-^100 = 4 X 10~^ ph cm~^ s~^and a photon index of 2.2. 



-46- 




Fig. 18. — Detection efficiencies in tlie LAT |b| >10° survey as a function of flux. Tfie solid line is 
for the entire blazar population while the dashed and long-dashed are for the FSRQs and BL Lacs 
respectively. The errors on the detection efficiency are due to the counting statistics in our Monte 
Carlo simulations. 
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Fig. 19. — Source count distribution for the whole extragalactic population (excluding the pulsars). 
The dashed line is the best power- law fit to the F(>100 MeV)> 7 x 10"^ ph cm s data. The 
inset shows the differential distribution. 
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Fig. 20. — Source count distribution for FSRQs. The dashed hne is the best power-law fit to the 
F(>100 MeV)> 7 x 10~^ph cm~^ data. The inset shows the differential distribution. 
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Fig. 21. — Source count distribution for BL Lacs. The dashed hne is the best power-law fit to the 
F(>100 MeV)> 3 x 10"^ ph cm s data. The inset shows the differential distribution. 




Fig. 22. — Luminosity functions of FSRQs in bins of redshift. The cumulative and differential 
distributions are shown, respectively, on the left and on the right panel. The (color-coded) solid 
lines are the ML fits to the 3 different datasets using a simple power law to model the GLF. 
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Fig. 23. — Luminosity functions of BL Lacs in bins of redshift. The cumulative and differential 
distributions are shown, respectively, on the left and on the right panel. The (color-coded) solid 
lines are the ML fits to the 2 different datasets using a simple power law to model the GLF. 



Table 1. The high-confidence association Bright AGN List 



FoM gtsrcid 



LAT Name 


Sourcs Name 


FoM 


Prob. 


SourcG Name 


Prob. 


Other Nam.es 






Class 


OFGL 


J0017. 


.4—0503 


CGRaBS J0017-0612 


16.20 





.93 


CGRaBS J0017-0512 





.92 




0, 


,227 


FSRQ 


OFGL 


J0033 


.6 — 1921 










BZB .70033-1921 





.99 


KUV 00311-1938 





610 


BLLac 


OFGL 


J0050 


.5 — 0928 


CGRaBS ,10050-0929 


61.01 


1 


.00 


CRATES J0050-0929 





.99 


PKS 0048-097 






BLLac 


OFGL 


J0051 


.1 — 0647 


CGRaBS ,70051-0650 


42.95 





.98 


CRATES J0051-0650 





.99 


PKS 0048-071 


1 


975 


FSR.Q 


OFGL 


J0112 


.1 + 2247 


CGRaBS ,70112+2244 


48.96 





.98 


S2 0109+22 


1 


.00 


S2 0109+22 





.265 


BLLac 


OFGL 


J0118. 


.7—2139 


CGRaBS J0118-2141 


35.02 





.97 


CGRaBS J0118-2141 





.99 


PKS 0116-219 


1 


,165 


FSRQ 


OFGL 


J0120 


.5—2703 


CGRaBS J0120-2701 


60.25 


1 


.00 


PKS 0118-272 


1 


.00 


PKS 0118-272 


0, 


.667 


BLLac 


OFGL 


J0136 


.6+3903 


BZB J0136+3905 


12.45 





.91 


B3 0133+388 


1 


.00 


B3 0133+388 






BLLac 


OFGL 


J0137. 


.1+4751 


CGRaBS J0136+4751 


52.21 





.99 


CGRaBS ,70136+4751 





.99 


DA 55 





,859 


FSRQ 


OFGL 


J0144 


.5+2709 


CRATES ,70144+2705 


30.93 





.96 


GRATES ,70144+2705 





.58 


TXS 0141+268 






BLLac 


OFGL 


JO 145 


.1 — 2728 


CGRaBS ,70145-2733 


37.41 





.97 


CGRaBS J0145-2733 





.96 


PKS 0142-278 


1 


.148 


FSRQ 


OFGL 


J0204 


.8 — 1704 


CGRaBS ,70204-1701 


55.22 





.99 


CGRaBS J0204-1701 





.96 


PKS 0202-17 


1 


,740 


FSRQ 


OFGL 


J0210. 


.8—5100 


CGRaBS J0210-5101 


69.87 


1 


.00 


PKS 0208-512 


1 


.00 


PKS 0208-512 


1 


,003 


FSRQ 


OFGL 


J0217, 


.8+0146 


CGRaBS J0217+0144 


52.67 





.99 


CGRaBS J0217+0144 


1 


.00 


PKS 0215+016 


1, 


.716 


FSRQ 


OFGL 


J0220 


.9+3607 


CGRaBS J0221+3556 


7.66 





.89 


CGRaBS J0221+3656 





.95 


B2 0218+36 


0, 


.944 


Unc*** 


OFGL 


J0222 


.6+4302 


BZB ,70222+4302 


23.18 





.95 


3C 66A 


1 


.00 


3C 66A 


0, 


,444 


BLLac 


OFGL 


J0229 


,S — 3(i40 


BZQ ,70229-3643 


29.20 





.96 


BZQJ0229-3643 





.94 


PKS 0227-369 


2 


,115 


FSRQ 


OFGL 


J0238 


.6+1636 


CGRaBS ,70238+1636 


60. ,54 


1 


.00 


CGRaBS ,10238+1636 


1 


.00 


AO 0235+164 





940 


BLLac 


OFGL 


,70245 


,6—4656 


CRATES ,70246-4651 


23.23 





.95 


CRATES ,70246-4651 





.54 


PKS 0244-470 






Unc*** 


OFGL 


J0303. 


.7—2410 


CRATES J0303-2407 


9.32 





.90 


PKS 0301-243 


1 


.00 


PKS 0301-243 


0, 


.260 


BLLac 


OFGL 


J0320 


.0+4131 


CGRaBS J0319+4130 


33.67 





.97 


0316+413 


1 


.00 


NGC 1275 


0, 


.018 


RG 


OFGL 


J0334 


.1—4006 


CGRaBS J0334-4008 


63.24 


1 


.00 


PKS 0332-403 


1 


.00 


PKS 0332-403 






BLLac 


OFGL 


J0349 


.8—2102 


CGRaBS J0349-2102 


47.40 





.98 


CGRaBS J0349-2102 





.99 


PKS 0347-211 


2 


,944 


FSRQ 


OFGL 


J0428 


.7—3755 


CGRaBS J0428-3756 


54.09 


Q 


.99 


CGRaBS J0428-3756 


1 


.00 


PKS 0426-380 


1 


,112 


BLLac 


OFGL 


J0449. 


.7—4348 


CRATES J0449-4350 


5.52 





.81 


PKS 0447-439 


1 


.00 


PKS 0447-439 


0, 


,205 


BLLac 


OFGL 


J0457. 


.1 — 2325 


CGRaBS J0457-2324 


35.74 





.97 


CGRaBS J0457-2324 


1 


.00 


PKS 0454-234 


1, 


,003 


FSRQ 


OFGL 


J0507, 


.9+6739 


BZB J0607+6737 


4.74 





.76 


lES 0502+676 


1 


.00 


lES 0502+675 


0, 


.416 


BLLac 


OFGL 


J0516. 


.2—6200 


CGRaBS J0616-6207 


12.04 





.91 


CGRaBS J0516-6207 





.94 


PKS 0516-621 






Unc*** 


OFGL 


J0531. 


.0+1331 


CGRaBS J0630+1331 


65.48 


1 


.00 


CRATES J0530+1331 


1 


.00 


PKS 0528+134 


2, 


.070 


FSRQ 


OFGL 


J0538 


.8—4403 


CRATES ,70538-4405 


53.80 





.99 


BZB,70538-4405 





.99 


PKS 0537-441 


0, 


.892 


BLLac 


OFGL 


J0654 


.3+4513 


CGRaBS ,70654+4514 


42.13 




.98 


CGRaBS ,106,54+4514 


1 


.00 


B3 0650+463 


0, 


.933 


FSRQ 


OFGL 


J0654 


.3+5042 


CGRaBS J0654+5042 


49.98 





.99 


CGRaBS J0654+5042 


1 


.00 








Unc*** 


OFGL 


J0700 


.0—6611 


CRATES J0700-6610 


33.82 





.97 


CRATES J0700-6610 





.64 


PKS 0700-661 






Unc*** 


OFGL 


J0712, 


.g_|_5034 


CGRaBS J0712+5033 


44.20 





.98 


CGRaBS J0712+5033 





.99 








BLLac 


OFGL 


J0714. 


.2+1934 


CLASS J0713+193B 


20.54 





.94 










0, 


.634 


FSRQ 


OFGL 


J0719 


.4+3302 


CRATES J0719+3307 


14.33 





.92 


BZUJ0719+3307 





.89 


TXS 0716+332 


0, 


.779 


FSRQ 


OFGL 


J0722 


.0+7120 


CGRaBS J0721+7120 


66.40 


1 


.00 


CRATES J0721+7120 


1 


.00 


S5 0716+71 


0, 


,310 


BLLac 


OFGL 


,70738 


.2+1738 


CGRaBS J0738+1742 


25.45 





.95 


PKS 0735+17 


1 


.00 


PKS 0735+178 


0, 


,424 


BLLac 


OFGL 


J0818. 


.3+4222 


CGRaBS J0818+4222 


61.26 


1 


.00 


OJ 425 


1 


.00 


OJ 426 


0, 


.630 


BLLac 


OFGL 


J0824, 


.9+5551 


CGRaBS J0824+BBB2 


57.80 


0, 


.99 


CGRaBS J0824+5662 


0, 


.98 


TXS 0820+560 


1, 


.417 


FSRQ 


OFGL 


J0855, 


.4+2009 


CGRaBS J0864+2006 


8.67 


0, 


.90 


OJ 287 


0, 


.99 


OJ 287 


0, 


.306 


BLLac 


OFGL 


J0921, 


.2+4437 


CGRaBS J0920+4441 


13.49 


0, 


.92 


CGRaBS J0920+4441 


0, 


.95 


RGB J0920+446 


2, 


.190 


FSRQ 


OFGL 


J0948, 


.3+0019 


CGRaBS J0948+0022 


18.64 


0, 


.93 


CGRaBS J0948+0022 


0, 


.94 


PMN J0948+0022 


0, 


.686 


FSRQ 


OFGL 


J0957 


.6+5522 


CRATES J0957+5522 


50.91 





.99 


BZQJ0957+5522 





.96 


4C +55.17 


0, 


,896 


FSRQ 


OFGL 


J1012 


.9+2435 


CRATES ,71012+2439 


13.63 





.92 










1 


,805 


FSRQ 


OFGL 


J1015, 


.2+4927 


CGRaBS J1015+4926 


18.06 


0, 


.93 


lES 1011+496 


1 


.00 


lES 1011+496 


0, 


.212 


BLLac 


OFGL 


J1015, 


.9+0515 


CRATES J1016+0513 


28.86 


0, 


.96 


CRATES J1016+0513 


0, 


.78 


PMN J1016+0612 


1, 


.713 


FSRQ 


OFGL 


J 1034, 


.0+6051 


CGRaBS J1033+60B1 


52.57 


0, 


.99 


CGRaBS J1033+6061 


0, 


.98 


S4 1030+61 


1, 


.401 


FSRQ 


OFGL 


J1053, 


.7+4926 


BZB J1063+4929 


11.55 


0, 


.91 


MS 1050.7+4946 


1 


.00 


MS 1060.7+4946 


0, 


.140 


BLLac 


OFGL 


J1054, 


.5+2212 


CLASS J1064+2210 


16.20 


0, 


.93 














BLLac 


OFGL 


J1057 


.8+0138 


CGRaBS J1058+0133 


3.71 





.68 


CGRaBS ,71058+0133 


0, 


.93 


PKS 1055+018 


0, 


,888 


FSRQ 


OFGL 


J 1058 


.9+5629 


CGRaBS ,71058+5628 


24.66 





.95 


R,XS ,110586+5628 


1 


.00 


RXS J105S6+5628 





,143 


BLLac 


OFGL 


JllOO, 


.2-8000 


CGRaBS J1058-8003 


53.65 


0, 


.99 


CGRaBS J1058-8003 


0, 


.99 


PKS 1057-79 






BLLac 


OFGL 


J1104, 


.5+3811 


CGRaBS J1104+3812 


35.10 


0, 


.97 


Mrk 421 


1 


.00 


Mrk 421 


0, 


.030 


BLLac 


OFGL 


J1129, 


.8-1443 


CRATES J1130-1449 


27.54 


0, 


.96 


BZQ J1130-1449 


0, 


.84 


PKS 1127-14 


1, 


,184 


FSRQ 


OFGL 


J1146 


.7-3808 


CGRaBS ,71147-3812 


45.04 





.98 


CGRaBS J1147-3812 





,99 


PKS 1144-379 


1 


,048 


FSRQ 


OFGL 


,11 159 


.2+2912 


CGRaBS ,71159+2914 


39.38 





.97 


CGRaBS J1159+2914 





.98 


4C 29.45 





.729 


FSRQ 


OFGL 


,11218 


.0+3006 


CGRaBS ,71217+3007 


31.80 





.96 


B2 1215+30 


1 


.00 


B2 1215+30 





,130 


BLLac 


OFGL 


,11221 


.7+2814 


CGRaBS ,71221+2813 


36.82 





.97 


W Com 


1 


.00 


W Com 





.102 


BLLac 


OFGL 


J1229, 


.1+0202 


CGRaBS J1229+0203 


73.53 


1 


.00 


3C 273 


1 


.00 


3C 273 


0, 


,158 


FSRQ 


OFGL 


J 1246, 


.6-2544 


CGRaBS J1246-2647 


43.45 


0, 


.98 


CGRaBS J1246-2547 


0, 


.99 


PKS 1244-255 


0, 


.636 


FSRQ 


OFGL 


J1253, 


.4+5300 


CRATES J1263+B301 


43.34 


0, 


.98 


S4 1260+53 


1 


.00 


S4 1260+53 






BLLac 


OFGL 


J1256, 


.1-0547 


CGRaBS J1266-0647 


71.21 


1 


.00 


3C279 


1 


.00 


3C 279 


0, 


.636 


FSRQ 


OFGL 


J1310, 


.6+3220 


CGRaBS J1310+3220 


55.91 


0, 


.99 


CGRaBS J1310+3220 


0, 


.99 


B2 1308+32 


0, 


.997 


FSRQ 
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Table 1 — Continued 



FoM gtsrcid 



LAT Name 


bource Name 


FoM 


Prob. 


Source Name 


Prob. 


Otiier Names 






(-^lass 




J1325 


.4 — 4303 


BZU J1325-4301 


75. 


23 




,00 


NGC 5128 




,00 


NGC 5128, Cen A 


0. 


,002 


RG 




J1331 


.7 — 0506 


CGRaBS J1332-0509 


of 


64 





,98 


CGRaBS .71332-0509 




.93 


PKS 1329-049 


2 


,150 


FSRQ 








CLASS ,11333+5057 






I' 


,94 










1 


,362 


FSR,Q 




J1355 




CRATES J1354-1041 


■ 


52 





,94 


BZUJ1354-1041 





.84 


PKS 1352-104 


0. 


,330 


FSRQ 




J1427 


.1+2347 


CRATES J1427+2347 


19. 


.69 





,94 


PKS 1424+240 




.00 


PKS 1424+240 






BLLac 




J1457 


.6 — 3538 


CGRaBS J1457-3539 


26. 


.03 





,95 


CGRaBS J1457-3539 





.99 


PKS 1454-354 


1. 


,424 


FSRQ 




J1504 


.4+1030 


CGRaBS J1504+1029 


48. 


.85 





,98 


CGRaBS J1504+1029 


^ 


.00 


PKS 1502+106 


1. 


.839 


FSRQ 


OFGL 


J1511 


.2 — 0536 


PKS 1608-05 


10. 


27 





,90 


BZQJ1510-0543 





.73 


PKS 1508-05 


1, 


.185 


FSRQ 


OFGL 


.11512 


.7-0905 


PKS 1510-08 


74, 


,49 


1 


,00 


BZQJ1512-0905 





.98 


PKS 1510-08 


0, 


.360 


FSRQ 


OFGL 


.T1517 


.9-2423 


CGR.aBS ,11517-2422 


19. 


18 





,94 


AP Lib 


1 


.00 


AP Lib 


0, 


.048 


BLLac 


OFGL 


J1522, 


.2+3143 


CGRaBS J1522+3144 


51. 


,06 


0, 


,99 


CGRaBS J1522+3144 


1 


.00 


TXS 1520+319 


1, 


.487 


FSRQ 


OFGL 


J1543, 


.1+6130 


CRATES J1542+6129 


45. 


,22 


0, 


,98 


RXS J15429+6129 


1 


.00 


RXS J15429+6129 






BLLac 


OFGL 


J1553, 


.4+1255 


CRATES J1553+1256 


26. 


,38 


0, 


,95 


PKS 1551+130 


0, 


.85 


PKS 1551+130 


1, 


,308 


FSRQ 


OFGL 


.T1555 


.8+1110 


CGRaBS ,71555+1111 


44, 


23 





,98 


PG 1553+11 


1 


.00 


PG 1553+11 





360 


BLLac 


OFGL 


.T1625 


.8-2527 


CGRaBS ,71625-2527 


56. 


82 





,99 


PKS 1622-253 





.99 


PKS 1622-253 





786 


FSRQ 


OFGL 


.J1635 


.2+3809 


CGRaBS ,71635+3808 


.54, 


10 





,99 


CRATESJ1635+3808 





.99 


4C +38,41 


1 


814 


FSRQ 


OFGL 


J1653 


.9+3946 


CGRaBS ,71653+3945 


59. 


,08 





,99 


Mrk 501 


1 


.00 


Mrli 501 





(133 


BLLac 


OFGL 


J1719, 


.3+1746 


CGRaBS J1719+1745 


40. 


,87 


0, 


,98 


PKS 1717+177 


1 


.00 


PKS 1717+177 


0, 


,137 


BLLac 


OFGL 


J17S1, 


.5+0935 


CGRaBS J1751+0939 


19. 


,73 


0, 


,94 


CGRaBS J1751+0939 


0, 


.99 


OT 081 


0, 


.322 


BLLac 


OFGL 


J1802, 


.2+7827 


CGRaBS J1800+7828 


28. 


,07 


0, 


,96 


CGRaBS J1800+7828 


0, 


.99 


S5 1803+78 


0, 


,680 


BLLac 


OFGL 


.T1847 


.8+3223 


CGRaBS ,71848+3219 


12, 


,76 





,92 


CGRaBS .71848+3219 





.94 


TXS 1846+322 


0, 


,798 


FSRQ 


OFGL 


.T1849 


.4+6706 


CGRaBS ,71849+6705 


53. 


,89 





,99 


CGRaBS .71849+6705 


1 


.00 


S4 1849+67 


0, 


,657 


FSRQ 


OFGL 


.J1911 


.2-2011 


CGRaBS ,71911-2006 


23. 


,51 





,95 


CGRaBS .71911-2006 





.97 


PKS 1908-201 


1 


,119 


FSRQ 


OFGL 


J1923 


.3-2101 


CGRaBS ,71923-2104 


37, 


72 





,97 


CGRaBS .71923-2104 





.97 


TXS 1920-211 





,874 


FSRQ 


OFGL 


J2000, 


.2+6506 


CGRaBS J1959+6508 


19. 


12 


0, 


,94 


lES 1959+650 


1 


.00 


lES 1959+650 


0, 


,047 


BLLac 


OFGL 


J2009, 


.4-4850 


CGRaBS J2009-4849 


72. 


,13 


1, 


,00 


PKS 2005-489 


1 


.00 


PKS 2005-489 


0, 


,071 


BLLac 


OFGL 


J2025, 


.6-0736 


CRATES J2025-0735 


42. 


,71 


0, 


,98 


BZQJ2025-0735 


0, 


.98 


PKS 2022-07 


1, 


,388 


FSRQ 


OFGL 


.T2056 


.1-4715 


CGRaBS ,72056-4714 


67, 


,00 


1 


,00 


CRATES ,72055-4716 


1 


.00 


PKS 2052-47 


1 


,491 


FSRQ 


OFGL 


.T2139 


.4-4238 


CRATES ,72139-4235 


13. 


,48 





,92 


MH 2136-428 


1 


.00 


MH 2136-428 






BLLac 


OFGL 


.T2143 


.2+1741 


CGRaBS ,72143+1743 


36. 


,88 





,97 


CGRaBS ,72143+1743 





.96 


OX 169 


0, 


.213 


FSRQ 


OFGL 


.T2147 


.1+0931 


CGRaBS ,72147+0929 


53. 


,97 





,99 


CGRaBS ,72147+0929 





.99 


PKS 2144+092 


1, 


.113 


FSRQ 


OFGL 


J2157, 


.5+3125 


CGRaBS J2157+3127 


54. 


,48 


0, 


,99 


CGRaBS J2157+3127 


0, 


.97 


B2 2155+31 


1, 


.486 


FSRQ 


OFGL 


J2158, 


.8-3014 


CGRaBS J2158-3013 


54. 


,87 


0, 


,99 


CGRaBS J2158-3013 


1 


.00 


PKS 2155-304 


0, 


.116 


BLLac 


OFGL 


J2202, 


.4+4217 


BZB J2202+4216 


45. 


,62 


0, 


,98 


BZB J2139-4239 


1 


.00 


BL Lacertae 


0, 


.069 


BLLac 


OFGL 


.12203 


.2+1731 


CGRaBS ,72203+1725 


23. 


,91 





,95 


CGRaBS ,72203+1725 





.93 


PKS 2201+171 


1, 


,076 


FSRQ 


OFGL 


.T2207 


.0-5347 


CGRaBS ,72207-5346 


39. 


56 





,97 


CGRaBS ,72207-5346 





.99 


PKS 2204-54 


1 


,215 


FSRQ 


OFGL 


.T2229 


.8-0829 


CGRaBS ,72229-0832 


42, 


,99 





.98 


CGRaBS J2229-0832 





.99 


PHL 5225 


1 


,560 


FSRQ 


OFGL 


.T2232 


.4+1141 


BZQ .72232+1143 


45, 


97 





.98 


BZQ ,72232+1143 


1 


.00 


CTA 102 


1 


,037 


FSRQ 


OFGL 


J2254, 


.0+1609 


CGRaBS J2253+1608 


70. 


,34 


1 


.00 


CGRaBS J2253+1608 


1 


.00 


3C 454.3 


0, 


,859 


FSRQ 


OFGL 


J232B, 


.3+3959 


CRATES J2325+3957 


29. 


,25 


0, 


.96 


B3 2322+396 


1 


.00 


B3 2322+396 






BLLac 


OFGL 


J2327, 


.3+0947 


CGRaBS J2327+0940 


21. 


12 


0, 


.94 


CGRaBS J2327+0940 


0, 


.93 


PKS 2325+093 


1, 


,843 


FSRQ 


OFGL 


J234B, 


.5-1559 


CGRaBS J2345-1555 


30. 


19 


0, 


.96 


CGRaBS J2345-1555 


0, 


.93 


PMN J2345-1555 


0, 


.621 


FSRQ 



***A11 these source have a flat radio spectrum but there are no other data reported in literature which allow to classify them 
either as FSRQs or BL Lacs. 
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Table 2. The low-confidence association Bright AGN List 



LAT Namr 



OFGL J0100.2+07B0 
OFGL J0238.4+2S55 
OFGL .70407.6-3829 
OFGL ,70412.9-5341 
OFGL ,70423.1-0112 
OFGL J0909. 7+0145 
OFGL J1034. 0+6051 
OFGL J1248. 7+5811 
OFGL J1625. 9-2423 
OFGL Jlfi41.4+.'-!<).'-!<( 
0FC-;L ,72(117,2 I ()(i()2 



FoM 



SoTirro Namr 



CRATES J0100+0746 
CGRaBS J0237+2848 
CRATES ,70406-3826 
CRATES ,70413-5332 
CGRaBS J0423-0120 
CGRaBS J0909+0200 
CRATES J1032+6051 

CRATES J1627-2426 

CLASS J1641+.'-!<(3.-> 
CL.VSS .I21Ji7 I ()(i():i 



FoM Proh 



6.12 0.78 

7.67 0.89 

3.00 0.61 

1.92 0.46 

4.26 0.72 

4.16 0.71 

5.22 0.79 

2.33 0.53 

6.22 t).nr, 

7,n:i n.f^s 



gtsrcid 



SoTirrr Namr 



CRATES J0423-0120 
PKS 0907+022 



Proti. Othrr Namrs 



Class 



0.000 

CGRaBS J0237+2848 0.88 B2 0234+28 1.213 

PKS 0405-385 1.285 

0.00 

0.84 PKS 0420-014 0.915 

0.87 PKS 0907+022 • ■ ■ 



PG 1246+586 



0.86 



1.064 



Unc"* 
FSRQ 

Unc*** 

Unc*** 

FSRQ 

BLLac 

FSRQ 

BLLac 

Unc*** 

FSRCJ 

7^iu'"'* 



***A11 these sources have a flat radio spectrum but there are no other data reported in literature which allow to classify 
them either as FSRQs or BL Lacs. 
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Table 3. The Fermi -LAT \b\ > 10° Bright AGN List 



LAT Name R.A. Dec ( 6 \/TS r» J^ioo'' Ppeak' ^25"^ Var. 



OFGL 


J0017 


4-0503 


4 


358 


5 


054 


101 


273 


—66 


485 


14 


7 


2 


71+0 


14 


13 


9 + 2 


4 


34 


8 + 6 


5 


12.1 + 1.4 


T 




J0033 


6-1921 


8 




— 19 


360 


94 




— 81 




10 


7 




70 i 


if 


J; 


6 i 






y ± i 


3 


0.4 i 0.1^ 


■_ 




J0050 


5-0928 




637 


— 9 


470 




209 


— 72 


341 


20 


5 




15 i 


08 


10 


-1- 1 
J It i 




19 


it 4 





8.8 it 1.3 






J0051 


1-0647 


io 


796 


-6 


794 




751 


-69 


666 


15 


^ 




22 + 






5 + 1 




19 


7 + 4 


^ 


7.2 + L4 




OFGL 


.101 00 


2+0750 


15 


051 


7 


844 


126 


716 


— 54 


963 


11 


1 


1 


8t) + 


16 


1 


9 + 


7 


3 


9+1 


7 


0.3 ± 0.1^ 




OFGL 


J0112 


1 + 2247 


18 


034 


22 


789 


129 


148 


— 39 


832 


17 


6 


2 


10 ± 


09 


7 


4 + 1 


2 


12 


6 + 2 


7 


6.0 + 0.7 




OFGL 


J0118 


7-2139 


19 


676 


— 21 


656 


172 


990 


— 81 


728 


17 


8 


2 


32 ± 


10 


9 


6 + 1 


4 


21 


4 + 4 


5 


7.6 + 1.1 


T 


OFGL 


J0120 


5-2703 


20 


128 


— 27 


056 


213 


951 


—83 


529 


11 


8 


1 


99 + 


14 


3 


2 + 


8 


6 


7 + 2 


3 


2.6 ± 0.8 




OFGL 


J0136 


6-1-3903 


24 


163 


39 


066 


132 


446 


— 22 


969 


12 


5 


1 


65 ± 


13 


1 


8 it 


5 


4 


7± 1 


5 


0.5 ±0.1''" 




OFGL 


J0137 


1-1-4751 


24 


285 


47 


854 


130 


818 


— 14 


317 


18 


8 


2 


20 ± 


09 


10 


9 it 1 


7 


18 


6 it 4 


5 


10.8 it 1.6 


T 


OFGL 


,70144 


5 + 2709 


26 


142 


27 


159 


137 


248 


— 34 


231 


10 


4 


2 


22 i 


14 


5 


4+1 


3 


12 


7 + 3 


8 


2.0 + 0.5 








1-2728 




289 










— 78 




"'"'^ 


^ 




55 ± 


14 




2+1 




26 


3 + 5 




9.4 + 1.3 


T 














068 


186 


072 


— 70 


274 




() 










1 + 1 






9 + 3 


9 


10.7 +1.3 








s'^ioo 










?Io 




~Va 


IIq 








28 + 






4 + 2 




7^ 








T 




J0217 


0-1-0140 


-iA 
















21 






13 + 






2 + 1 






5 + 3 


a 




T 


OFGL 


J 0220 


y-|-oDU ( 


35 


243 


36 


121 


142 


504 


23 


325 


12 






61 +; 


16 




i 2 


4 


22 


5 it 6 




10.9 +1.3 




OFGL 


J0222 




35 


653 


43 


043 


140 


132 





763 


47 






97 it 


04 


25 


9 it 1 


Q 


49 


6 i 4 


g 


26 6 i 1 4 


X 


OFGL 


J0229 


5-3640 


37 


375 


— 36 


681 


243 


801 


fjY 


189 


10 


2 


2 


57 + 




15 


8 + 2 




34 


1 + 6 


2 


14 1 + 1 5 


T 


OFGL 


J 023 8 


Ji_-j-2855 


39 


600 


28 


923 


149 


521 


28 


368 


10 




2 


49 i 


15 


9 


+ 2 


() 


24 


7 + 5 


9 


8.6 + 1.6 






































6 + 2 










6 7 . () + 2.2 


T 


OFGL 


7024^ 




41 


423 


46 


934 


262 


019 


~60 


098 


11 




^ 


34 + 


1^ 


"^6 


2 + 1 




^12 


4 + 4 








OFGL 


J uouo 


7-2410 


45 


940 


— 24 


176 


214 


764 


QQ 


119 


12 


3 


2 


01 i 


13 


3 


8 + 


g 




+ 2 


g 


2 9 + 09 




OFGL 




0+4131 


50 


000 


41 


524 


150 


601 


13 


230 


29 




2 


17 i 


06 


22 


1 + 1 


9 


35 


9 it 5 


3 


18.2 + 1.4 


T 


OFGL 


J0334 


1-4006 


53 


546 


4Q 


107 


244 


710 


— 54 


088 


13 


2 


2 


15 + 


12 


5 


3 + 1 


1 


11 


2 + 3 


1 


4.9 + 1.4 




OFGL 


J0349 


8-2102 


57 


465 


— 21 


046 


214 


385 


4g 


035 


21 


2 


2 


55 i 


09 


19 


2 + 2 


3 


27 


8 + 5 





17 3 + 16 




OFGL 


J 04.07 


6-3829 


01 


923 


— 38 


491 


241 


360 


— 47 


751 


13 




2 


_|_ 


13 


7 


5 + 1 


5 


22 


2 + 4 


\ 


6.9 + 1.3 


T 


OFGL 


JO412 




63 


230 


— 53 


686 


263 




44 


716 


10 












4+1 


3 


12 


3 + 3 




6 + 13 




OFGL 


J 04 23 


1-0112 


65 


785 


_ \ 


204 


195 


131 


— 33 


092 


\ \ 




2 


38 ± 


16 




1 + 2 


2 


13 


4 + 4 





10.5 + 3.1 




OFGL 


J0428 


7-3755 


67 


193 


— 37 


923 


240 


689 


—43 


597 


39 


Q 


2 


14 + 


05 


24 


5 + 1 


g 


31 


5 + 4 


7 


23.1 ± 1.6 




OFGL 


J 0449 


7-4348 


72 


435 


43 


815 


248 


780 


—39 


859 


28 


4 


2 


01+0 


06 


12 


0+1 


3 


21 


1 + 4 


2 


12.2 + 1.4 




OFGL 


J 045 7 


1-2325 


74 


288 


23 


432 


223 


739 


34 


880 


52 


3 


2 


23 i 


04 


41 


8 i 2 


3 


64 


2 + 6 


4 


OD.D It J- .0 


X 


OFGL 


J0507 


y-|-D ( oy 


76 


985 


67 


650 


143 


772 


15 


905 


13 




^ 


67 i 


18 




7 i 


® 




-1- 1 


^ 


0.3±0.lt 




OFGL 


J0516 


conn 


79 


063 


— 62 


000 


271 


376 


— 34 


834 








17 i 


17 


5 


A 4- ^ 


7 




1 -1- Q 


3 


0.4 lb 0. 1 ' 




OFGL 


J0531 


0+1331 


82 


761 


13 


528 


191 


385 


— 10 


992 


17 


3 


2 


54 ± 


09 


24 


3 + 2 


9 


39 


5 + 6 


7 


00 is _i_ 1 

ZO.O ± Z. I 


T 


OFGL 
















ool 








6 


^ 








1' ^ ^ 


^ 


= 1^ 


7 + 5 








OFGL 


J0654 


3+4513 


103 


590 


45 


220 


171 


228 


19 


369 


29 


2 




32 i 


06 


23 


8 + 2 




56 


4 + 7 





20.3 i 1.6 


■p 


OFGL 


J0654 


3+5042 


103 


592 


50 


711 


165 


676 


21 


107 


15 


6 


2 


00 ± 


10 


5 


5 it 1 


1 


9 


5 it 2 


6 


4.9 ±1.3 


X 


OFGL 


J0700 


0-6611 


105 


016 


—66 


199 


276 


778 


—23 


809 


10 


1 


1 


98 ±: 


14 


3 


9 ± 1 





8 


7 it 2 


6 


0.4 ± 0.1 ^ 




OFGL 


J0712 


9+5034 


108 


231 


50 


575 


166 


688 


23 


900 


11 


2 


2 


04 + 


14 


3 


9± 1 


1 


10 


5 it 2 


7 


3.3 it 0.7 




OFGL 


J0714 


2+1934 


108 


552 


19 


574 


197 


685 


13 


648 


15 





2 


35 it 


10 


10 


7 it 1 


6 


27 


it 5 





10.0 it 1.6 


X 


OFGL 


J0719 


4+3302 


109 


869 


33 


037 


185 


139 


19 


855 


12 


3 


2 


37 + 


15 


7 


8 + 1 


7 


20 


8 + 4 


9 


7.5 it 1.5 


X 


OFGL 


JO 722 


0+7120 


110 


508 


71 


348 


143 


976 


28 


029 


34 


4 


2 


08 ± 


05 


16 


4+1 


4 


29 


+ 4 


2 


17.0 it 1.6 


X 


OFGL 


J0738 


2+1738 


114 


575 


17 


634 


201 


933 


18 


081 


11 


9 


2 


10 + 


14 


4 


6 + 1 


1 


7 


5 + 2 


4 


3.6 it 1.4 






J0818 






579 




367 


178 




33 


409 


20 


9 




07 it 


08 


9 


It i 


3 


14 


5±2 


9 


7.0 it 1.1 




OFGL 


J0824 


9+5551 


126 


239 


55 


859 


161 


981 


35 


142 


10 


6 


2 


81 ±0 


20 


11 


4 + 2 


9 


42 


OitS 


1 


10.8 ± 1.3 


X 


OFGL 


J0855 


4+2009 


133 


857 


20 


162 


206 


810 


35 


974 


15 


1 


2 


31 ±0 


11 


9 


Oil 


5 


19 


0it4 


1 


7.8 it 1.3 




OFGL 


J0909 


7+0145 


137 


446 


1 


757 


228 


640 


31 


262 


11 


6 


2 


67it0 


16 


10 


4 + 2 


1 


22 


9it6 


1 


9.5 it 0.3 




OFGL 


J0921 


2+4437 


140 


320 


44 


617 


175 


809 


44 


876 


15 


2 


2 


35 + 


12 


8 


6 + 1 


5 


15 


7 + 4 


2 


9.2 it 1.4 




OFGL 


J0948 


3+0019 


147 


077 





317 


236 


530 


38 


549 


12 


8 


2 


60 + 


14 


12 


1 + 2 


2 


29 


2 + 5 


7 


9.1 it 1.4 


X 


OFGL 


J0957 


6+5522 


149 


424 


55 


375 


158 


605 


47 


939 


24 





2 


01+0 


07 


8 


7+1 


1 


12 


9 + 3 





9.2 it 1.3 




OFGL 


J1012 


9+2435 


153 


241 


24 


598 


207 


897 


54 


406 


12 


4 


2 


22 + 


12 


6 


lit 1 


2 


10 


9it3 


6 


4.5 it 0.9 


X 


OFGL 


J1015 


2+4927 


153 


809 


49 


463 


165 


473 


52 


727 


23 


8 


1 


73it0 


07 


4 


9it0 


7 


7 


lit 1 


7 


8.9it 1.6 




OFGL 


J1015 


9+0515 


153 


991 


5 


254 


236 


457 


47 


036 


20 


6 


2 


20 ±0 


08 


11 


7it 1 


5 


21 


8±4 


5 


13.1 it 1.5 


X 


OFGL 


J 1034 


0+6051 


158 


504 


60 


853 


147 


765 


49 


122 


14 


8 


2 


48it0 


13 


9 


3it 1 


7 


22 


0it4 


7 


7.5 it 1.3 




OFGL 


J1053 


7+4926 


163 


442 


49 


449 


160 


309 


58 


263 


10 


1 


1 


42 + 


20 





5 + 


3 


1 


8 + 


9 


0.2±0.lt 




OFGL 


J1054 


5+2212 


163 


626 


22 


215 


216 


968 


63 


049 


11 


2 


2 


24 + 


15 


4 


9 + 1 


3 


10 


8 + 3 


3 


4.3 it 1.0 




OFGL 


J1057 


8+0138 


164 


451 


1 


643 


251 


219 


52 


709 


10 


3 


2 


20 + 


17 


5 


+ 1 


4 


10 


7 + 2 


8 


9.2 ±1.7 




OFGL 


J1058 


9+5629 


164 


731 


56 


488 


149 


521 


54 


442 


12 





2 


11±0 


14 


3 


9±1 





8 


3±2 


7 


5.0± 1.8 




OFGL 


JllOO 


2-8000 


165 


057 


-80 


012 


298 


047 


-18 


212 


12 


1 


2 


71 ±0 


16 


17 


1±3 


8 


38 


4 + 8 


5 


11.1±2.2 


X 


OFGL 


J 1104 


5+3811 


166 


137 


38 


187 


179 


868 


65 


056 


47 


1 


1 


77 ±0 


04 


15 


3±1 


1 


20 


9it3 


1 


15.9 ± 1.3 




OFGL 


J1129 


8-1443 


172 


454 


-14 


727 


275 


133 


43 


694 


10 


5 


2 


69it0 


18 


9 


9it2 


4 


25 


8it 5 


8 


10.8 ±1.6 




OFGL 


J1146 


7-3808 


176 


689 


-38 


149 


289 


170 


22 


988 


10 


4 


2 


21+0 


14 


5 


7+1 


4 


7 


5 + 2 


8 


3.5±1.2 




OFGL 


J1159 


2+2912 


179 


800 


29 


216 


199 


605 


78 


307 


14 


6 


2 


47 + 


13 


10 


3 + 1 


8 


16 


+ 3 


8 


9.4 ± 1.0 




OFGL 


J1218 


0+3006 


184 


517 


30 


108 


188 


826 


82 


097 


27 


4 


1 


89 + 


06 


9 


7+1 


1 


40 


9 + 4 


7 


10.4 ± 1.0 


X 


OFGL 


J1221 


7+2814 


185 


439 


28 


243 


201 


593 


83 


336 


24 





1 


93it0 


07 


8 


3it 1 


1 


17 


2 + 3 


5 


7.5 ±0.9 


X 


OFGL 


J1229 


1+0202 


187 


287 


2 


045 


289 


975 


64 


355 


52 





2 


71it0 


05 


75 


2 + 4 


3 


137.0 ± 13 





65.5 ±2.6 


X 



- 55 - 



Table 3 — Continued 



LAT Name 


R.A. 


Dec I 


b 




pa 


-r'lOO 






Var. 



OFGL 


J1246, 


,6-2544 


191 


,655 


-25, 


,734 


301, 


.571 


37. 


125 


11 


7 


2 


.24 ±0. 


,14 


6. 


.8 + 1. 


.6 


15.3 + 4, 


.3 


7. 


.6 + 1. 


,4 




OFGL 


J1248. 


7+5811 


192 


,189 


58, 


,191 


123. 


.617 


58. 


.934 


14. 


.3 


1 


.95 ±0. 


,11 


3. 


.8 + 0. 


.8 


8.0 + 2, 


.4 


7. 


.4+1. 


,6 




OFGL 


J1253, 


,4+5300 


193, 


,369 


53, 


,001 


122, 


229 


64. 


125 


12 


.1 


2 


.17±0. 


,14 


4. 


.7+1. 


.1 


9.1 + 2, 


,6 


5. 


.6 + 1. 


,5 




OFGL 


J1256, 


,1-0547 


194, 


,034 


-5, 


,800 


305, 


.081 


57. 


.052 


36. 


,8 


2 


.35 ±0. 


,05 


31. 


.5 + 2. 


,3 


46.3 + 6, 


,8 


29. 


.7+1. 


,8 


T 


OFGL 


J1310, 


,6+3220 


197, 


,656 


32, 


,339 


85, 


.458 


83. 


.331 


27. 


,3 


2 


.25 ±0. 


,07 


15. 


.5 + 1. 


,6 


37.3 + 4, 


,6 


16. 


.4+1 


,1 


T 


OFGL 


J1325, 


,4-4303 


201 


,353 


-43, 


,062 


309, 


.501 


19. 


.376 


12 


,4 


2 


.91 ±0. 


,18 


21. 


.5 + 4. 


,5 


32.3 + 8, 


,0 


22. 


2 + 2 


,4 




OFGL 


J1331 


,7-0506 


202 


,935 


-5, 


,112 


321, 


247 


56. 


.320 


14. 


,3 


2 


.59 ±0. 


,12 


13. 


.0 + 2. 


,1 


33.0 + 5, 


,9 


10. 


.7+1 


2 


T 


OFGL 


J1333, 


,3+5058 


203, 


,331 


50, 


,973 


107, 


.300 


64. 


.865 


12 


4 


2 


.40 ±0. 


,14 


7. 


.2 + 1. 


,5 


13.7 + 4, 


,6 


9. 


1 + 1 


,3 


T 


OFGL 


J1355, 


,0-1044 


208, 


,764 


-10, 


,735 


327, 


221 


49. 


113 


11 


,5 


2 


.37 ±0. 


,15 


7. 


.6 + 1. 


,8 


34.4 + 5, 


,5 


8. 


.7+1 


,3 


T 


OFGL 


J1427, 


,1+2347 


216 


,794 


23, 


,785 


29, 


.472 


68. 


166 


24. 


,1 


1 


.80 ±0. 


,07 


6. 


.2 + 0. 


,8 


8.7 + 2, 


,2 


5. 


1 + 1 


,0 




OFGL 


J1457, 


,6-3538 


224, 


,407 


-35, 


,639 


329, 


.936 


20. 


.530 


39. 


,6 


2 


.24 ±0. 


,05 


36. 


.6 + 2. 


,4 


77.2 + 7, 


,1 


32. 


1 + 0. 


,5 


T 


OFGL 


J1504, 


,4+1030 


226 


,115 


10, 


,505 


11. 


.409 


54. 


577 


88. 


,2 


2 


.17±0. 


,02 


81. 


.4 + 2 


,7 


260.0 + 15, 


.0 


69. 


.3 + 2. 


,1 


T 


OFGL 


J1511 


,2-0536 


227, 


,814 


-5, 


,613 


354. 


.099 


42. 


.948 


10. 


,8 


2, 


.41 + 0. 


,15 


8. 


.8 + 2. 


,1 


16.2 + 4, 


.4 


8. 


.6+1 


,7 




OFGL 


J1512 


,7-0905 


228 


,196 


-9, 


,093 


351. 


282 


40. 


153 


45. 


,0 


2, 


.48 + 0. 


,05 


55. 


.8 + 3. 


,3 


165.9+ 11, 


,7 


50. 


.6 + 2. 


,3 


T 


OFGL 


.11517, 


,9-2423 


229 


,496 


-24, 


,395 


340. 


724 


27. 


521 


12 


,3 


1 


.94 + 0. 


,14 


4. 


1 + 1. 


,2 


7.0 + 2, 


,4 


5. 


.2 + 0. 


,6 




OFGL 


J1522 


,2+3143 


230, 


,552 


31, 


,726 


50. 


143 


57. 


014 


34. 


,3 


2, 


.39 + 0. 


,06 


25. 


7 + 2 


,1 


42.0 + 5, 


,1 


22. 


.2+1 


,5 


T 


OFGL 


J1543, 


,1+6130 


235, 


,784 


61, 


,504 


95. 


.383 


45. 


.370 


10. 


,5 


2 


.00 + 0. 


,15 


2. 


.5 + 0. 


,7 


4.3 + 1, 


,7 


3. 


.3+1 


,4 




OFGL 


J1553, 


,4+1255 


238, 


,368 


12 


,922 


23. 


746 


45. 


225 


23. 


,7 


2 


.23 + 0. 


,07 


16. 


1 + 1 


,8 


33.6 + 5, 


,6 


15. 


.6 + 2. 


,2 


T 


OFGL 


J1555, 


,8+1110 


238, 


,951 


11 


,181 


21. 


.911 


43. 


941 


31. 


,5 


1 


.70 + 0. 


,06 


8. 


.0+1. 


,0 


11.6 + 2, 


,3 


10. 


2 + 2 


,0 




OFGL 


J1625, 


,8-2527 


246 


,470 


-25, 


,451 


352. 


164 


16. 


.308 


11 


,4 


2 


.40 + 0. 


,15 


16. 


.0 + 4. 


,6 


28.4 + 8, 


,0 


19. 


.8+1 


,3 




OFGL 


J 1625. 


9-2423 


246 


,494 


-24, 


,393 


353. 


.005 


16. 


.995 


10. 


,1 


2 


.46 + 0. 


,14 


19. 


.9 + 5. 


,1 


32.1 + 8, 


,1 


10. 


.7 + 0. 


,9 




OFGL 


J1635, 


,2+3809 


248 


,821 


38, 


,158 


61. 


118 


42. 


.333 


27. 


,3 


2 


.44 + 0. 


,07 


22. 


.0 + 2. 


,2 


49.8 + 6, 


,0 


17. 


.6+1 


,3 


T 


OFGL 


J 1641. 


4+3939 


250, 


,355 


39, 


,666 


63. 


.239 


41. 


.239 


17. 


,7 


2 


.43 + 0. 


,10 


13. 


.4 + 2 


,1 


33.7 + 6, 


,3 


12. 


.7+1 


,4 


T 


OFGL 


J1653, 


,9+3946 


253, 


,492 


39, 


,767 


63. 


612 


38. 


841 


19. 


,0 


1 


.70 + 0. 


,09 


3. 


1 + 0. 


,6 


6.9 + 1 


,8 


3. 


.3 + 0. 


,8 






J1719 


,3+1746 


259 


,830 


17 


,768 


39. 


553 


28. 


080 


23 


,3 




.84 i 


,07 


6. 


9 + 


,9 


15.0 + 3 


,0 


6. 


.3+1 


,2 




OFGL 


J1751 


,5+0935 


267, 


,893 


9, 


,591 


34. 


.867 


17. 


614 


23. 


,1 


2, 


.27 + 0. 


,07 


18. 


.4 + 2 


,1 


41.4 + 6 


,3 


17. 


.8+1 


,9 


T 


OFGL 


J1802, 


,2+7827 


270, 


,567 


78, 


,466 


110. 


.026 


28. 


.990 


12 


,6 


2, 


.25 + 0. 


,14 


6. 


.0 + 1. 


,4 


11.1 + 3, 


,1 


5. 


.9+1 


,3 




OFGL 


J1847, 


,8+3223 


281 


,954 


32, 


,385 


62. 


.065 


14. 


.838 


16. 


,0 


2, 


.37 + 0. 


,10 


14. 


.7 + 2. 


,4 


28.0 + 4, 


,9 


9. 


.4 + 0. 


,6 


T 


OFGL 


J1849, 


,4+6706 


282 


,365 


67, 


,102 


97. 


.503 


25. 


027 


28. 


,0 


2, 


.17 + 0. 


,06 


15. 


.9 + 1. 


,5 


28.8 + 4, 


,1 


14. 


.9+1 


,5 


T 


OFGL 


J1911 


,2-2011 


287, 


,813 


-20, 


,186 


16. 


.818 


-13. 


.266 


20. 


,0 


2, 


.43 + 0. 


,08 


22. 


.5 + 2. 


7 


52.3 + 7, 


,2 


18. 


.7 + 0. 


,8 


T 


OFGL 


J1923, 


,3-2101 


290, 


,840 


-21 


,031 


17. 


.205 


-16. 


199 


16. 


,4 


2 


.31 + 0. 


10 


13. 


1 + 2. 


,0 


41.6 + 6, 


,1 


11. 


.3 + 0. 


,6 


T 


OFGL 


J2000, 


,2+6506 


300, 


,053 


65, 


,105 


97. 


974 


17. 


.630 


15. 


,8 


1 


.86 + 0. 


11 


4. 


.2+1 


,0 


6.3 + 2, 


,1 


3. 


.4+1 


,3 




OFGL 


J2009, 


,4-4850 


302, 


,363 


-48, 


,843 


350. 


.361 


-32. 


.607 


10. 


,9 


1 


.85 + 0. 


,14 


2. 


.9 + 0. 


,9 


5.5 + 2, 


,1 


3. 


.0 + 0. 


,6 




OFGL 


J201 7. 


2+0602 


304, 


,302 


6, 


,048 


48. 


.596 


-15. 


.991 


12 


,7 


1 


.87 + 0. 


,12 


3. 


.7 + 0. 


,9 


6.6 + 2, 


,3 


0.6 + 0.1+ 




OFGL 


J2025 


,6-0736 


306, 


,415 


-7, 


,611 


36. 


.883 


-24. 


.389 


50. 


,6 


2, 


.30 + 0. 


,04 


48. 


.0 + 2. 


,6 


73.6 + 7, 


,1 


43. 


.0 + 2. 


,0 


T 


OFGL 


J2056, 


,1-4715 


314, 


,034 


-47, 


,251 


352. 


.586 


-40. 


.358 


12 


,5 


2, 


.56 + 0. 


,15 


11. 


1 + 2. 


,3 


21.1 + 5, 


,2 


10. 


.7+1 


,7 




OFGL 


J2139, 


,4-4238 


324, 


,865 


-42 


,642 


358. 


237 


-48. 


.332 


20. 


,1 


2, 


.01 + 0. 


,08 


8. 


.0 + 1. 


,2 


13.1 + 3, 


,0 


7. 


.7+1 


,3 




OFGL 


J2143 


,2+1741 


325, 


,807 


17, 


,688 


72. 


.016 


-26. 


.051 


14. 


5 


2 


.57 + 0. 


12 


14. 


1 + 2. 


,2 


30.7 + 6, 


,2 


12. 


.0+1 


7 




OFGL 


J2147, 


,1+0931 


326, 


,777 


9, 


,519 


65. 


.805 


-32. 


.236 


19. 


,9 


2 


.53 + 


10 


16. 


.6 + 2. 


1 


34.7 + 6, 


,1 


16. 


.6+1 


,6 


T 


OFGL 


J2157, 


,5+3125 


329, 


,384 


31, 


,431 


84. 


747 


-18. 


.258 


10. 


,0 


2 


.41 + 0. 


,15 


7. 


.5+1. 


,7 


13.9 + 3, 


,9 


7. 


.3+1 


,5 




OFGL 


J2158, 


,8-3014 


329, 


,704 


-30, 


,237 


17. 


711 


-52. 


.236 


43. 


,9 


1 


.85 + 0. 


,04 


18. 


1 + 1. 


,2 


29.2 + 3, 


,6 


18. 


.5+1 


,4 


T 


OFGL 


J2202 


,4+4217 


330, 


,622 


42 


,299 


92. 


.569 


-10. 


.398 


12 


,3 


2, 


.24 + 0. 


,12 


8. 


.5 + 1. 


,8 


12.8 + 4, 


,3 


8. 


.0 + 2. 


,0 




OFGL 


J2203 


,2+1731 


330, 


,815 


17, 


,532 


75. 


.715 


-29. 


.529 


12 


,7 


2, 


.25 + 0. 


,13 


6. 


.9 + 1. 


,4 


17.0 + 3, 


,5 


8. 


.3+1 


,6 


T 


OFGL 


J2207, 


,0-5347 


331 


,765 


-53, 


,786 


339. 


.948 


-49. 


.832 


12 


,4 


2, 


.65 + 0. 


,17 


11. 


.5 + 2. 


,5 


54.6 + 8, 


,0 


11. 


1 + 1 


,8 


T 


OFGL 


J2229 


,8-0829 


337, 


,452 


-8, 


,495 


55. 


.326 


-51. 


701 


16. 


,8 


2 


.67 + 0. 


12 


15. 


.9 + 2. 


,4 


27.7+5 


,7 


12. 


.0 + 0. 


,4 




OFGL 


J2232 


,4+1141 


338, 


,117 


11 


,690 


77. 


372 


-38. 


.592 


15. 


,2 


2 


.61 + 0. 


,12 


14. 


.0 + 2. 


,3 


24.6 + 6, 


,2 


11. 


.2+1 


,3 




OFGL 


J2254, 


,0+1609 


343, 


,502 


16, 


,151 


86. 


125 


-38. 


187 


149 


,1 


2, 


.41 + 0. 


,02 


246. 


1 + 5. 


,2 


385.8 + 20, 


,5 


221. 


6 + 4. 


,3 


T 


OFGL 


J2325 


,3+3959 


351 


,334 


39, 


,993 


105. 


.532 


-19. 


.952 


11 


,4 


1 


.89 + 0. 


,13 


2. 


.8 + 0. 


,8 


11.0 + 2, 


,7 


1. 


.3 + 0. 


,4 


T 


OFGL 


J2327, 


,3+0947 


351 


,833 


9, 


,794 


91. 


159 


-47. 


821 


17. 


,1 


2, 


.73 + 0. 


,12 


18. 


.3 + 2. 


,6 


51.0 + 8, 


,4 


15. 


.8+1 


,6 


T 


OFGL 


J2345, 


,5-1559 


356, 


,389 


-15, 


,985 


65. 


677 


-71. 


.092 


15. 


,5 


2, 


.42 + 0. 


,12 


10. 


.5 + 1. 


,7 


22.3 + 4, 


,3 


10. 


.3+1 


,3 


T 



^Spectral index derived from a single power-law fit over the 0.2-100 GeV energy range 

'^Flux {E > 100 MeV, in 10~^ ph cm~^ s""'") derived from a single power-law fit over the 0.2-100 GeV energy range 
''Weekly averaged peak flux (E > 100 MeV) in 10~^ ph cm~^ s""^ 

'^Flux {E > 100 MeV, in 10~^ ph cm~^ ^~^) obtained by adding the fluxes estimated in the two energy ranges 0.1— 1 GeV and 1— 100 GeV 

^ Flux at E > 1 GeV in 10 ph cm"'^ s""*" . For these sources, only an upper limit is obtained for the 0.1 — 1 GeV flux fsee |Abdo et all i2009d )). 
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Table 4. Sources in both Fermi-LAT and EGRET samples 



LAi JName 


tj(jKhl IName 


mean 


peak 


index 


100 


-pE 
^ peak 


index 


Type 


OFGL J0210.8— 5100 


J0210— 5055 


24.3 


76.2 


2.28 


85.5 


134.0 


1.99 


FSRQ 


OFGL J0222. 6+4302 


J0222+4253 


25.8 


49.6 


1.96 


18.7 


25.3 


2.01 


BLLac 


OFGL J0238.6+1636 


J0237+1635 


72.5 


104. 


2.05 


25.9 


65.1 


1.85 


BLLac 


OFGL J0423. 1—0112 


J0422— 0102 


8.0 


13.3 


2.37 


16.3 


81.7 


2.44 


FSRQ 


OFGL J0457. 1—2325 


70456—2338 


41.7 


64.2 


2.23 


8.1 


18.8 


3.14 


FSRO 


OFGL J0516. 2-6200 


J0512-6150 


5.3 


11.1 


2.17 


7.2 


28.8 


2.40 


Unc 


OFGL J0531 0+1331 


J0530+1323 


24.3 


39.4 


2.54 


93.5 


351.0 


2.46 


FSRQ 


OFGL J0538.8— 4403 


J0540-4402 


37.6 


49.6 


2.18 


25.3 


91.1 


2.41 


BLLac 


OFGL J0722. 0+7120 


J0721+7120 


16.3 


29.0 


2.07 


17.8 


31.8 


2.19 


BLLac 


OFGL J0738. 2+1738 


J0737+1721 


4.6 


7.46 


2.10 


16.4 


37.5 


2.60 


BLLac 


OFGL J0855. 4+2009 


J0853+1941 


8.9 


18.9 


2.30 


10.6 


15.8 


2.03 


BLLac 


OFGL J092 1.2+4437 


J0917+4427 


8.6 


15.6 


2.34 


13.8 


40.8 


2.19 


FSRQ 


OFGL J0957.6+5522 


J0952+5501 


8.7 


12.8 


2.01 


9.1 


47.2 


2.12 


FSRQ 


OFGL J1104.5+3811 


J1104+3809 


15.3 


20.9 


1.76 


13.9 


27.1 


1.57 


BLLac 


OFGL J1159. 2+2912 


J1200+2847 


10.3 


16.0 


2.47 


7.5 


163.0 


1.73 


FSRQ 


OFGL J1229. 1+0202 


J1229+0210 


75.2 


136. 


2.71 


15.4 


53.6 


2.58 


FSRQ 


OFGL J1256. 1—0548 


71255—0549 


31.5 


46.3 


2.34 


74.2 


267.0 


1.96 


FSRO 


OFGL J1325.4— 4303 


J1324-4314 


21.4 


32.3 


2.90 


13.6 


38.4 


2.58 


RG (CenA) 


OFGL J1333. 3+5058 


J1337+5029 


7.2 


13.7 


2.4 


9.2 


26.8 


1.83 


FSRQ 


OFGL J1457. 6-3538 


J1500-3509 


36.5 


77.2 


2.24 


10.9 


40.7 


2.99 


FSRQ 


OFGL J1512.7— 0905 


J1512-0849 


55.8 


165. 


2.47 


18.0 


51.1 


2.47 


FSRQ 


OFGL J1517.9-2423 


J1517-2538 


4.1 


6.96 


1.93 


8.4 


53.3 


2.66 


BLLac 


OFGL J1625. 8-2527 


J1626-2519 


16.0 


28.4 


2.39 


21.3 


90.2 


2.21 


FSRQ 


OFGL J1625. 9-2423 


J1627-2419 


19.9 


32.1 


2.45 


15.8 


55.2 


2.21 


Unc 


OFGL J1635.2+3809 


J1635+3813 


22.0 


49.7 


2.43 


58.4 


108.0 


2.15 


FSRQ 


OFGL J1802.6-3939'' 


J1800-3955 


25.4 


64.0 


2.23 


9.8 


189.0 


3.10 


Unc 


OFGL J1833.4-2106'' 


J1832-2110 


42.0 


56.8 


2.61 


26.6 


99.3 


2.59 


FSRQ 


OFGL J1911. 2-2011 


J1911-2000 


22.4 


52.2 


2.42 


17.5 


47.6 


2.39 


FSRQ 


OFGL J1923. 3-2101 


J1921-2015 


13.0 


41.6 


2.31 


4.6 


31.0 


2.10 


FSRQ 


OFGL J2025.6-0736 


J2025-0744 


47.9 


73.5 


2.30 


21.2 


74.5 


2.38 


BLLac 


OFGL J2056.1-4715 


J2055-4716 


11.0 


21.0 


2.55 


9.6 


35.0 


2.04 


FSRQ 


OFGL J2158. 8-3014 


J2158-3023 


18.0 


29.1 


1.85 


13.2 


30.4 


2.35 


BLLac 


OFGL J2202.4+4217 


J2202+4217 


8.4 


12.8 


2.23 


11.1 


39.9 


2.60 


BLLac 


OFGL J2232.4+1141 


J2232+1147 


14.0 


24.5 


2.61 


19.2 


51.6 


2.45 


FSRQ 


OFGL J2254.0+1609 


J2254+1601 


246. 


385. 


2.41 


53.7 


116.0 


2.21 


FSRQ 



■"lO-^ph cnr^s^i 
''source located at |&| < 10° 
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Table 5. Correlation analysis for the radio and gamma-ray properties 



Data 


Method Source type 


# objects 


Correlation coeff. 


Chance probability'^ 


<S'8.4GHz — -F'> 100 MeV, peak 


Spearman rank All 


106 


0.42 


0.0000045 


Log Lr-hog i-y 


Sp. r., partial All 


90 


0.46 


0.0000025 


'S'8.4GHz — -F>100MeV,peak 


Spearman rank FSRQ 


57 


0.19 


0.080 


Log Lr-Log 


Sp. r., partial FSRQ 


57 


0.34 


0.0047 


-58.4 GHz — -F'> 100 MeV, peak 


Spearman rank BL Lacs 


42 


0.49 


0.00055 


Log Lr-Log 


Sp. r., partial BL Lacs 


30 


0.60 


0.00023 



"^if the samples were unbiased 
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Table 6. Composition of the |b| >10° sample (TS>100). 

CLASS # objects 

Total 132 

FSRQs 57^ 

BL Lacs 42^ 

Uncertain'^ 5^ 

Radio Galaxies 2^ 
Pulsars 

Anti-associations 4 

Low confidence associations 10 

Unassociated sources 5 

'^Part of the high confidence sample. 

^Blazars with uncertain classification. 

''Five LAT detected pulsars plus OFGL 
J0025.1-7202 (47 Tuc) and OFGL J0538.4- 
6856 (asso ciated with t he Lar ge Magellan 



Cloud, see lAbdo et al.l (|2009d ^). 
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Table 7. Results of the best power-law fits to the source counts distributions. Errors within 
brackets are systematic uncertainties due to the incompleteness of the sample. The lower part of 
the table shows the results for the flux-limited portion of the sample. 



SAMPLE 


# Objects 


a 




EDB fraction"^ 


AlF 


121 


2.59±0.12 


2.62±0.24 


7.2% 


Blazars 


106 


2.50±0.12 


2.24±0.22(±0.24) 


6.1% 


FSRQs 


57 


2.60±0.14 


2.15±0.28(±0.32) 


3.1% 


BL Lacs 


42 


2.34±0.15 


0.41±0.06(±0.06) 


1.0% 


FSRQs 


29 


2.52±0.20 


1.93±0.35(±0.09) 


2.6 %^ 


BL Lacs 


9 


2.50±0.37 


0.48±0.16(±0.02) 


1.3% 



'^In units of 10^ cm^ s deg ^. 

'^Fraction of the EGRET diffuse extragalactic background 
( Sreekumar et al.lll998l ) resolved into sources by LAT for 4 x 10~^ <Fioo < 
10~^ph cm~^ s~^. 

'^Includes all sources except 7 pulsars and 4 anti-associated objects. 

^The lower limit of integration in Eq. [7] has been set to 6 x 10~^ ph cm~^ 
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Table 8. Results of the Y/Ymax test. 



SAMPLE 


# Objects 


<V/Vmax> 


FSRQs 


57 


0.645 ±0.043 


BL Lacs 


31 


0.430 ±0.055 


BL Lacs^ 


42 


0.472 ±0.046 


BL Lacs^ 


42 


0.473 ±0.046 


All with z>0 


92 


0.512 ±0.031 


FSRQs " 


29 


0.654 ±0.061 


BL Lacs ^ 


8 


0.542 ±0.103 



'^For all the 11 BL Lacs without redshift, we 
have assumed z = <z> where <z>=0.38. 

^For all the 11 BL Lacs without redshift, 
we have drawn a redshift measurement from 
a Gaussian distribution with mean 0.38 and 
dispersion 0.34. 

'^Flux-limited sample. 

Tabic 9. Results of best-fit power-law models to the GLFs of FSRQs in different redshift bins. 
The lower part of the table shows the results for the flux-limited portion of the sample. Errors 
represent the 68 % confidence level. Uncertainties within bracketes are systematic errors due to 

the incompleteness of the sample. 



Redshift bin 


# Objects 


f3 


Normalization^ 


z = 0.0 - 


0.9 


20 


1.57±0.10 


2.43±0.52(±0.36) 


z = 0.9 - 


1.4 


17 


2.56±0.29 


5.59±1.33(±0.83) 


z > 1.4 




20 


2.58±0.19 


13.07±2.92(±1.96) 


z = 0.0 - 


0.9 


10 


1.46±0.13 


1.73±0.54(±0.08) 


z = 0.9 - 


1.4 


8 


2.65±0.48 


6.62±2.34(±0.33) 


z > 1.4 




11 


2.63±0.30 


17.76±5.35(±0.88) 



^Normalization of the GLF model at 10^^ erg s ^ expressed 
in units of 10~^^ erg~^ sMpc~^. 
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Table 10. Results of best-fit power-law models to GLFs of BL Lacs in dilferent redshift bins. 
The lower part of the table shows the results for the flux-limited portion of the sample. Errors 
represent the 68 % confidence level. Uncertainties within bracketes are systematic errors due to 

the incompleteness of the sample. 



Redshift bin 


# Objects 


(3 


Normalization^ 


z = 0.0 - 0.3 


15 


2.08±0.16 


7.67±1.98(±3.06) 


z > 0.3 


16 


2.10±0.11 


4.75±1.18(±1.90) 


z > 0.0 


31 


2.17±0.05 


4.19±0.75(±1.67) 


z > 0.0 


8 


1.90±0.11 


6.70±2.36(±0.33) 



'^Normalization of the GLF model at 10^* erg s ^ expressed 
in units of 10~^^ erg~^ sMpc~^. 



